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GENETIC MEASUREMENT OF THE HALF LIFE OF Bi?! 


J. SosNIAK? AND R. E. BELL 


ABSTRACT 


The half life of Bi?°? has been measured by the genetic method, using the 
parent—daughter pair Po?*? and Bi?°7. Both activities are electron capturers, and 
the measured half life depends only on the activity ratio between them. The 
value found, 71/2 = 28+3 years, agrees with one of three previous determinations 
of this half life. 


INTRODUCTION 


It is necessary to know the half life of Bi?°? in order to make measurements 
of activation cross sections of reactions leading to Bi?®’, such as Pb? (p, 27). 
Previous measurements of this half life by different workers had led to dis- 
cordant values. We now report the result of a measurement by a genetic 
method, leading to the result 7), = 28+3 years. 

In the first detailed work on Bi?°7, Neumann and Perlman (1951) estimated 
the half life both >y a genetic method involving the a-decay of At?", and by 
Geiger counting sources of Bi?®? over a period of 2.75 years. The first method 
vielded a value of 50 years, and the second, showing no perceptible decay, 
suggested a lower limit of 40 years. Neumann and Perlman therefore quoted 
the value 50 years without assigning a probable error. 

Harbottle (1955 and private communication) used a balanced double ioniza- 
tion chamber, standardized with radium, to observe the decay of Bi?’ over 
a period of 3 months. His value for the half life was 27+(2 or 3) years. 

Cheng et al. (1955) observed the decay of a sample of long lived bismuth 
activity for a period of 11 years, and found a half life of 8.0+0.6 years. A 
scintillation spectrum of this source appeared to be ‘identical with that of 
the 50 year Bi?®’”’. This result clearly disagrees with the results of the other 
investigators. 


EXPERIMENTAL METHOD 
The present measurement emerges from a study of the radioactive chain 
Po?*’—EC—Bi?*?—EC-—Pb?” (stable). 
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A measured amount of pure Po?®’ is allowed to decay completely (71,2. = 5.7 
hours), and the amount of Bi?” activity that grows from it is measured. The 
ratio of the two activities is equal to the inverse ratio of their half lives, and 
the half life of Bi?®’ is thus determined. The essence of the experiment is just 
the measurement of the ratio of two similar electron capture activities. In the 
actual experiment the K X-rays emitted by the two activities were counted 
with a sodium iodide crystal and a 28-channel analyzer. The K X-rays form 
a very clean peak (actually two peaks because of the escape effect in Nal) 
in the scintillation spectrum, easily separable from the effects of any gamma 
rays. 

The intrinsic efficiency of the sodium iodide counter is close to unity for 
the K X-rays, and it only remains to estimate the number of K X-rays emitted 
per disintegration of Po?®’ and Bi®®’, This question has been discussed for a 
number of electron capturers, including Po?®’, by Bell and Skarsgard (1956). 
Their conclusion was that the ratio of K capture to total electron capture 
from all shells could not be far from 1:1.33. This conclusion could be upset 
only for an electron capturer whose decay energy was low enough to be 
comparable with the A electron binding energy of the atom, but in the case 
of Po?’ an estimate exists for the total decay energy, 2.80+0.45 Mev (Gray 
1955). We thus conclude that there are 0.75 K shell vacancies per electron 
capture process in Po**?, Additional K shell vacancies may arise from A 
conversion of gamma rays following the electron capture process. If a com- 
plete table of gamma ray intensities and conversion coefficients were at hand, 
this contribution could be worked out. Alternatively, as shown by Bell and 
Skarsgard, the contribution may be measured by detecting K—K X-ray coin- 
cidences from the source by detecting the ‘‘sum peak” in the scintillation 
spectrum of a weak source placed close to the counter. Their measurement 
for Po®®? showed that the number of K X-rays was increased by 9% owing 
to internal conversion. The total number of K shell vacancies per disintegration 
of Po?’ therefore becomes 0.82. 

For Bi?®’, a complete disintegration scheme is available from the work of 
Alburger and Sunyar (1955), including a close estimate of the decay energy ; 
for one of the electron capture branches this energy is so low that only L 
and higher shell capture occurs. With the help of the compilation of Brysk 
and Rose (1955), the ratio of A to total electron capture is found to be 
1: 1.54 for the sum of all the electron capture branches, leading to the figure 
of 0.65 K vacancies per disintegration. The effect of K conversion of gamma 
rays is computed from Alburger and Sunyar’s work to add 0.11 to this figure, 
giving a total of 0.76 A vacancies per disintegration. The counting efficiencies 
for Bi? and Po?’ are therefore very nearly equal, being in the ratio 0.76:0.82. 
Effects of Auger transitions cancel between the two activities, and absorption 
within the sources is eliminated from consideration by making the sources of 
Bi?®? and Po?” physically identical. 

The Po?’ parent material was produced by the reaction Bi? (p, 3m) in the 
McGill synchrocyclotron. Bismuth trioxide powder was wrapped in 0.0005 in. 
aluminum foil to form a strip-like target 0.007 in. thick, and bombarded in 
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the circulating beam with about 1 microampere-hour of 24 Mev protons. 
According to the cross-section curves of Bell and Skarsgard (1956), the 
maximum yield would have occurred at about 30 Mev, but the lower energy 
was preferred in order to minimize the formation of 9-day Po? by the (p, 4n) 
reaction. The bombardment also produces some Bi?” directly by the (p, p2n) 
reaction; the chemical procedure described below removed this component. 
In addition, there would be appreciable production of the 2.9-year Po?®’ by 
the (p, 2m) reaction; fortunately, however, Po?’ emits relatively few X-rays 
and its presence does not affect our measurements. 

A few minutes after the end of the bombardment, the Po*®? was separated 
from the active bismuth trioxide powder by a simple chemical procedure. The 
objects of this separation were to get rid of the bulk of the bombarded material, 
and to avoid the Bi?” formed during the bombardment either by the (p, p2z) 
reaction or by the decay of Po*®’. The active powder was dissolved in hydro- 
chloric acid to which a few milligrams of tellurium carrier had been added. 
The addition of a few drops of stannous chloride precipitated the tellurium, 
carrying the Po?” activity with it. The Po?’ separation was considered to 
have taken place at the moment of formation of this precipitate. The preci- 
pitate was then centrifuged off and redissolved in hydrochloric acid. This 
simple procedure sufficed for the purpose because most of the detectable 
active atoms present in the target were Po?” in the first place. 

On account of the large ratio of the half lives, about 40,000:1, the initial 
Po?*? sample was much too strong to be counted in direct comparison with 
the eventual Bi?’ daughter. Instead, the Po? solution was diluted three times 
by successive factors of 10, producing a nominal dilution factor of 1000. 
Polonium activity is notorious for its tendency to adhere to glass surfaces, and 
in the dilution procedure polyethylene utensils were used, except for glass 
graduated pipettes. Each step in the dilution was checked by counting the 
Po?’ activities, and the counting results, rather than the volumetric measure- 
ments, were always accepted as representing the true dilution. The discrepancy 
was sometimes as high as 16%. The result of the dilution procedure was a 
set of four identical polyethylene bottles each holding 5 ml of solution, the 
first one containing 0.9 of the original amount of Po?®’, and the last (nomin- 
ally) 10-* of the original amount. The actual counting comparison then 
involved the Po?’ in the weakest sample and the Bi?’ formed by the decay 
of the strongest sample; the activity ratio between these two samples was 
about 40:1, so that they could be counted in the same geometry. 

The weakest Po?*? sample was counted at once in direct contact with a 
1.5in.X1lin. sodium iodide crystal. Its 5.7-hour decay was followed for a 
few hours, and the counting rate was extrapolated back to the moment of 
chemical separation. The strong source was allowed to stand for several days 
to allow the Po?’ to decay to Bi?®’, and then counted in the same manner as 
the Po?°?, The Bi?®? sample was checked for the absence of the 9-day Po? 
over a period of some days, and checked again after some weeks or months 
to show that the Bi?*? sample was not showing any decay faster than that 


-expected for a half life of several years. 
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RESULTS 

Three complete experiments were performed giving three numerical values 
for the half life of Bi®®’, viz. 24.3, 29.7, and 29.2 years. The average of the 
three values is 27.7 years with a root mean square deviation in the average 
of +1.4 years. This deviation is of little value as an estimate of the real error 
in the result. We have considered the errors in the various factors entering 
into the experiment, and estimate a probable error of about 10% for them 
all. This error would be higher if it were not for the way in which many 
uncertain factors cancel for the two sources being compared. We therefore 
give for the result of the experiment, 


Ti (Bi?®?) = 28+3 years. 


This result agrees very well with that of Harbottle (1955), and is not in 
strong contradiction with the early work of Neumann and Perlman (1951). 
The value of Cheng et al. for the half life, 8 years, remains as before, 
unexplained. 
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POISSON BRACKETS IN FIELD THEORY! 


HANS FREISTADT 


ABSTRACT 


Poisson brackets for covariant field theory are defined in such a way as to 
demonstrate the close connection and ready transition between the classical 
brackets and the corresponding commutators of quantum theory. The approach 
of Good is followed in general; but the questions of tensor algebra are handled 
differently, requiring the introduction of a family of space-like surfaces and 
their normals. As an illustration, the free Klein-Gordon and Dirac fields are 
worked out. 


INTRODUCTION 


It is our purpose in this paper to develop a Poisson bracket (P.B.) formalism 
for classical covariant field theory, in such a way as to permit ready passage 
to quantum theory. Hamiltonian formalisms for classical field theory have 
been much neglected, except for the work of Bergmann (1949) and his dis- 
ciples, and that of Good (1954). Treatises on quantum field theory (e.g. 
Wentzel 1943; Schweber, Bethe, DeHoffmann 1955) rarely go beyond the 
Lagrangian formalism; then they recall that in quantum mechanics, com- 
mutation relations between momenta and co-ordinates are introduced, and 
proceed to do likewise for the fields. The P.B., the classical source of the 
commutation relations, are usually not even mentioned, let alone developed. 
More recent approaches (Jauch and Rohrlich 1955) forsake the Hamiltonian 
formalism altogether, deriving commutation relations from an action principle. 
However, if one wishes to write out explicitly the (functional) Schrédinger 
equation for the quantized field system, from there to pass to the Hamilton- 
Jacobi formulation and to the non-linear modification of the original field 
equations (Freistadt 1956), the development of a covariant Hamiltonian 
formalism is rather necessary. The recovery of the field equations from the 
equations of motion in quantum mechanical P.B. form is usually the only 
method available for ascertaining whether a proposed Hamiltonian density 
deserves confidence as a description of the system under consideration. This 
is especially so in those cases in which the Hamiltonian does not follow from 
the Lagrangian in the elementary way, as in the Dirac field. 

Bergmann (1949) does treat P.B. in detail; but his approach is based on 
a ‘parameter formalism”’, the effect of which is that field momentum densities 
are of the same tensor rank as the field variables to which they are conju- 
gate. Without implying any criticism of Bergmann’s approach, we do not 
follow it here; rather, we let our momentum densities be of one higher tensor 
rank than the corresponding field variables (four momentum symbols for 
every field variable), following in this respect the procedure first introduced 
by Born (1934) and Weyl (1984), and subsequently used by Schwinger (for 


1Manuscript received August 20, 1958. 
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which see Jauch and Rohrlich 1955) and Good (1954). The latter’s purely 
classical Hamiltonian formalism is entirely adequate, and is, in fact, somewhat 
more elegant than ours. However, Good’s formalism is not a convenient 
starting point for Hamiltonian quantization, in which P.B. are replaced by 
commutators (or anticommutators). To begin with, Good’s P.B., unlike com- 
mutators, are not antisymmetric, i.e. (u,v), #4 —(v, u),. Another difficulty 
derives from the tensor rank of Good’s P.B., (uw, v),, which, unlike the com- 
mutator, is one higher than the product of the variables involved, because of 
differentiation with respect to the momentum densities. This is avoided here 
by a contraction with the normal unit vector to a family of space-like surfaces, 
assumed to be given. In our procedure, the P.B. (say, of field variables and 
their conjugate momenta) have the usual values of quantum mechanical 
commutators (except for trivial factors); this permits almost automatic 
quantization. 

Our definition of P.B. is thus dependent on the particular family of space- 
like surfaces chosen. However, measurements (at least in flat space-time— 
very little is known about measurements in curved space) always involve 
quantities (such as momentum, energy) which are the result of an integration 
over a space-like surface. This is especially true in the quantized theory, in 
which measurements pertain to a system of particles. In the course of such an 
integration, the explicit dependence on the family of space-like surfaces dis- 
appears. This is the usual state of affairs in quantum theory. Good’s formalism 
is not tied to such a family of space-like surfaces; that is perhaps the reason 
which renders it unsuitable for quantization. 

Our formalism is entirely covariant and couched in the language of geo- 
metrical objects in space-time; the latter is not assumed to be flat. However, 
the components of the metric tensor are not treated as field variables; our 
treatment therefore does not cover the theory of gravitation. 

An earlier attempt (Freistadt 1957) based on a tensor Hamiltonian density, 
appears to be less satisfactory than the present modification of Good’s approach. 


NOTATION AND REVIEW OF EARLIER WORK 

As in our earlier papers (1955, hereafter quoted as CFHJ, and 1956, quoted 
as QFHJ), the field variables are denoted by y*; they are functions of the 
space-time co-ordinates &, which are not required to be Cartesian. The metric 
2Zuy is assumed to be given. The dynamical behavior of the fields is obtained, 
as usual, from a variational principle, |W = 0, where W = (1 ofL dw, and dw 
is a scalar volume element of integration in space-time. The Lagrangian 
density L(y*, y¥,, &”) is a scalar function. The symbol y*, means Vv, y*, the 
covariant derivative of the field variables. There is assumed to be given a 
continuous family of non-intersecting, but not necessarily plane, space-like 
surfaces. The symbol do, is an element of integration on one such surface. The 
unit normal at each space-time point to the space-like surface passing through 
that point is indicated by nu*. The momentum densities are defined by 


ry" = OL/dy*,, and the Hamiltonian density by 
(1 ) HH (y, r) =L — y4yr4%, 
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provided the elimination of the y*, is possible. (Hereafter, as in eq. (1), the 
indices are omitted where y, 7, £ appear as arguments.) 
The Hamiltonian equations of motion are those of Good, 


‘a yt, = —dH/an,!, 
Vat s* = 0H/dy* 


(except for a trivial change in sign, traceable to the definition of the Hamil- 
tonian, and consistent with earlier papers). 
CLASSICAL POISSON BRACKETS 

The classical P.B. of two variables, u(y, 7) and v(y, 7), at the same 
space-time point, is defined by 
(3) (u,v) = n* [(du/dme*)(dv/dy°) — (du/ dy") (dv/ Ame*)]. 
For instance, 
(4) (ra", vy?) = n"b,°. 
The P.B. so defined has all the usual algebraic properties (see, for example, 
Dirac 1947); in particuiar, (u,v) = —(v, 4), a result which is useful when 
passing to quantum theory. 


In the definition of P.B. of quantities at different points on a space-like 
surface, a functional derivative is used. It is defined as 


(5) du(t) /dv(t’) = Ou/ dv 6(E—€’) 


where du/dv is the formal derivative, and 6(—£’) is a three-dimensional 
delta function, meant to be integrated over a space-like surface. Note that, 
as usual, the functional derivative does not have the same dimensions as the 
formal derivative. 

We now introduce 


af u(t) du(t’) éu(é) —v(é") ) maa “ 
t) v(t’)) = a = Te 1 : 
(u(&), v(é’) | n (, 2 By"(e"”) By" (E”) be (E) n (£’’) dag(’’) 


(6) = J (u,v) 5(E—#') 5(e’— 2") n*(E"") doa (”), 


For instance, 
(7) (wa"(E), y?(E')) = 64%54(E—2’) 
where 


b«(E—£’) = nd (E—2*’), 


Sf(e’) 8% (E—&) dow (#’) = (8), 


as defined in QFHJ and also in Jauch and Rohrlich (1955). 
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EQUATIONS OF MOTION 
Proceeding as in Good (1954), one verifies directly that, K(y, x) being 
a dynamical field variable, 
0K /dy4 = na(ra4%, K) 
n*(0K/dr4*) = —(y*, K) 
from which it follows (by letting K = H, and referring to eq. (2)) that the , 
Hamiltonian equations of motion can be written in the form 
nry4, oa (y4, H) 
Vata* = Na (wa?, ED). 


(8) 


| 
' 


(9) 


The first of these equations is formally a contraction with n* of the first of 
eqs. (2); but it is identical with it, in view of the completely arbitrary choice 
of the family of space-like surfaces. Equations (9), in the quantum mechanical 
form (eqs. (12) below) are used in quantum field theory to test whether a 
Hamiltonian is the correct one to use. Equations analogous to (9) and (12) 
are well known in the literature of quantum field theory, though nowhere, as 
far as we know, do they appear in this completely covariant form. 


QUANTIZATION 
Passage to quantum theory is now effected in the usual Hamiltonian way 
by defining the quantum mechanical P.B. 


(10) the(u, V)am = [U, 


where the brackets on the right-hand side mean, as usual, the commutator 
[ ]_ or the anticommutator [| ]|,, to be used depending on the spin of the 
field in question, except that in the equations of motion (9) the quantum 
mechanical P.B., as is well known, is always the commutator. The quantum 
mechanical P.B. are then assigned exactly the same values as the classical 
P.B., which gives, for instance, for the commutator (or anticommutator) 


(11) [waX(E), y2(E’)] = thed,%64(E—£’), 


that is, the expression used in Jauch and Rohrlich (1955) and in QFHJ. 

Then the test for the suitability of a Hamiltonian density in quantum 
field theory consists in verifying whether eq. (9) in quantum mechanical 
form, i.e. 


nvy4, = (ihc) [y*, A] 
Vata% = (thc)—n.[x,%, A] 


(12) 


reproduce the field equations of the theory in question. As is well known, 
in applying the usual algebraic identities to the quantum mechanical P.B., 
the order of the factors, which are now quantum mechanical operators, must 
be kept in mind. In the two examples below, each term in the Hamiltonian 
density // consists of two factors, and the right hand side of eq. (12) is ex- 


panded by the familiar identity [a, bc]. = [a, b]4c—Ofc, a],. 








FREISTADT: POISSON BRACKETS IN FIELD THEORY 9 


THE FREE KLEIN-GORDON FIELD 
From the usual Lagrangian density (see, for example, CFHJ), one obtains, 
by eq. (1), the Hamiltonian density 


(13) Hgoa = dmc *y+2k(hc)—'9*a,. 


From the first of eq. (12), applied to y and y*, one obtains the usual definition 
of the momentum densities (see CFH]J), while the second of eq. (12), applied 
to m and x*, yields the Klein-Gordon field equations for y and y*. (Note 
that the result is (VaV*-+x)y = 0, since we use the signature (—1, —1, —1, 1).) 


THE FREE DIRAC FIELD 


The Hamiltonian density here cannot be obtained from the Lagrangian 
density (see CFHJ) by the usual method; but various conjectures suggest 


(14) Hp = —Hinwyay. 


That this is, in fact, a suitable Hamiltonian density is seen from the applica- 
tion of (12). Only one non-vanishing momentum density is left, obtained 
from the Lagrangian formalism (CFHJ). The first of eq. (12) leads to the 
Dirac equation (i7*Va+«)~¥ = 0, while the second leads to its adjoint 
tVaWaasy¥* —kKWaay = 0. 
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CANADIAN STANDARD OF FREQUENCY! 


S. N. Kaura,? C. F. Pattenson,? M. M. THomson‘ 


ABSTRACT 


Over the past 3 years a frequency standard of very high precision has been 
installed in Canada. It is composed of equipment located in three different 
laboratories in Ottawa, Ontario, but separated by a few miles. Intercomparison 
of frequency between these laboratories, which is done by sending signals over 
telephone lines and related techniques, is briefly described. Results indicate 
frequency stability of about 2:10!° over short and long periods. Absolute fre- 
quency is determined from astronomical observations. International inter- 
comparison is carried out by phase measurement of standard frequency and by 
observations of time signals; some of the results are presented. 


INTRODUCTION 


The interest in accurate measurement of frequency has increased recently 
because of the advances in atomic and molecular physics, and in radio com- 
munication and navigation. In physical measurements the need for accurate 
frequency determination and the measurement of time intervals is of long 
standing. Most national laboratories, including our own, have had frequency 
standards for a long time. However, the increase in accuracy to about 1:10" 
is comparatively recent. In this paper we will attempt a brief description of 
the steps taken in Canada to achieve a highly precise and stable standard of 
frequency. 

Since the definition of time is so intimately connected with the definition 
of frequency the definition of one without the other is not practicable. The 
National Research Council and the Dominion Observatory in Ottawa, there- 
fore, work in close co-operation. All astronomical observations are carried 
out at the Dominion Observatory. Three laboratories, namely, the Division 
of Applied Physics and the Radio and Electrical Engineering Division of the 
National Research Council and the Division of Positional Astronomy of the 
Dominion Observatory jointly compose and maintain the Canadian Standard 
of Frequency. 

APPARATUS 

In this paper features common to the three laboratories will be described. 
The apparatus may be divided into three parts: (1) crystal oscillators, (2) 
intercomparison and monitoring equipment, and (3) astronomical observations. 
The three laboratories are situated in Ottawa at a distance of 4 to 6 miles 
from one another. Intercomparison among these laboratories is made by 
transmitting standard frequency signals over telephone wires. 

1Manuscript received September 9, 1958. 

Contribution from the Divisions of Applied Physics and Radio and Electrical Engineering, 
National Research Council, and the Dominion Observatory, Ottawa, Canada. 

Issued as N.R.C. No. 4990. 

*National Research Council, Division of Applied Physics, Ottawa, Canada. 


8National Research Council, Division of Radio and Electrical Engineering, Ottawa, Canada. 
‘Dominion Observatory, Division of Positional Astronomy, Ottawa, Canada. 


Can. J. Phys. Vol. 37 (1959) 


10 








KALRA ET AL.: CANADIAN STANDARD OF FREQUENCY 11 


The crystal oscillators used as refereince oscillators at 100 kc each consist of 
an Essen type quartz ring cut in the Z plane and mounted on silk suspensions 
(Essen 1951), an oven and a resistance bridge type of oven control (Booth and 
Laver 1942), bridge stabilized oscillator (Booth and Laver 1946), and the 
required power supplies. These oscillators are protected against electrical 
mains failure by quick-starting emergency mains power supplies. All the 
crystal units were made by the British Post Office Radio Research Laboratories 
and the oscillators obtained earlier were also made by them. In the oscillators 
obtained at a later date the circuitry was made by Messrs. Airmec Limited 
under licence. These units are mounted so that they operate free from shock 
and vibration and are maintained at a constant ambient temperature. Since 
they are used as reference standards their frequencies are never adjusted and 
on some of the units the frequency trimmer capacitor has been left out. The 
physical distribution of these oscillators is as follows: Applied Physics Division 
has four, Radio and Electrical Engineering Division has three, and the 
Dominion Observatory has two. In addition to these units each laboratory 
has a number of other oscillators of various manufacture which are used 
for other purposes. 

Each laboratory has its own monitoring equipment. Usually two types of 
records are kept. The frequencies of all the oscillators are intercompared 


FREQUENCY INTERCOMPARISON OF QUARTZ RING OSCILLATORS 


a x10+9 





Fic. 1. Spot frequency measurements. Oscillators AP, BP, and FP are located at the 
Division of Applied Physics, 5DO at the Dominion Observatory, and R2R at the Division of 
Radio and Electrical Engineering. Measurements are made at the Division of Applied Physics 
laboratories. The actual fractional difference of frequency expressed in parts per 10° {(Af/f) X 10° 
= —} is listed in the top left corner of each graph. The notation on the left gives the scale in 
each case. 
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daily using integrating periods of a few seconds each, i.e. spot frequency 
measurements. The crystal oscillators used as clocks are also monitored for 
net time gained or lost between different clocks for each successive 24-hour 
period. From this the 24-hour average frequency is also derived. Beat counters 
of commercial manufacture are normally used for integrating the number of 
beats between different sets of crystal oscillators. Accuracy of intercomparison 
of crystal clocks is about 5:10" when using a 24-hour integrating period. 
This can be increased to about 5:10” if the need arises. The accuracy for spot 
frequency measurements is usually better than 1:10" and this can also be 
increased if the need arises. Figure 1 shows typical results of spot frequency 
measurements. 

The reference frequency signal obtained from the quartz crystal oscillators 
at 100 kc is divided down to 1 c.p.s. by electronic circuitry. The lower fre- 
quencies serve as reference time intervals in addition to being reference 
frequencies. The signals at 1 c.p.s. serve as clocks and the astronomical obser- 
vations and other standard time measurements are made with reference to 
them. Frequency multipliers are used for obtaining higher reference fre- 
quencies. Figure 2 shows a photograph of the main part of the installation at 
the Division of Applied Physics. 

In certain types of equipment, e.g. the oven for the crystal, it is essential 
to maintain continuity of a-c. mains supply. Such units are protected against 
mains failure by automatic emergency mains supplies with switching and 
starting times of a fraction of a second. This helps in maintaining stability 
of operation and continuity of record. The crystal oscillators are protected 
against variation of loading by means of isolating amplifiers. Other pre- 
cautions are taken to ensure stability of operation and long operating life of 


components. 


ASTRONOMICAL OBSERVATIONS 


The unit of time, the second, formerly was defined as a fraction of the 
period of rotation of the earth, and as such suffered the same irregularities 
that were inherent in the rate of rotation of the earth. It is now related to 
the time, which is the independent variable in the definition of the motion of 
the earth about the sun, and is designated as Ephemeris Time. The second 
of time is defined (C.I.P.M. 1957) as the fraction 1/31 556 925.974 7 of the 
tropical year 1900 January 0 at 1200 hours Ephemeris Time (E.T.). The year 
referred to is the year for which this date is the middle. For the physical 
realization of this definition two types of observations are usually made. They 
are (i) the meridian passage of the stars which yields Universal Time (U.T.2); 
and (ii) the position of the moon with respect to the stars which yields AT 
(see equation 1). At the Dominion Observatory, stellar transits are observed 
with a Photographic Zenith Telescope (Thomson 1955) and the lunar obser- 
vations are made with a Markowitz Moon Camera (Markowitz 1954; Thom- 
son 1958). Figure 3 shows a photograph of this equipment. U.T.2 is a good 
approximation to a uniform time scale, but is affected by small and un- 
predictable variations in the period of rotation of the earth which amount to 
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a few parts in 10° over short periods of time (Essen et al. 1958). Observations 
of the moon are made in U.T.2, and the resulting difference between the 
observed and predicted positions of the moon is interpreted as the accumulated 
error in U.T.2 with respect to E.T. Hence the relationship 

(1) AT = E.T.-—U.T.2. 

The quantity A7’ is being determined from star observations made at many 
observatories combined with the results of about 20 moon cameras built for 
the International Geophysical Year and dispersed around the world. The 
result is expected to give a more accurate definition of E.T. than was available 
with former techniques. 

The difference between a second of U.T.2 and the second of E.T. is small 
and usually the second of U.T.2 is used in the definition of frequency. In the 
future it will be possible to use the second of E.T. for more precise definition 


of frequency. 


INTERLABORATORY INTERCOMPARISON OF FREQUENCY 

As stated earlier the three component parts of the Canadian Frequency 
Standard are separated by a few miles. A signal at the fundamental frequency 
of 100 kc from one of the Essen ring quartz crystal oscillators is transmitted 
to the Applied Physics Division laboratories over standard Bell Telephone 
lines from each of the other two laboratories. There are no repeater amplifiers 
on the way and the lines are not disturbed or otherwise used by the Bell 
Telephone Company. At the Applied Physics Division laboratories these 
signals are amplified and used in the same manner as the local standard fre- 
quency signals. Because of the length of line involved the spot frequency 
measurements have a slightly higher uncertainty, about 1:10° for the signal 
from the Dominion Observatory and about 5:10° for the signal from Radio 
and Electrical Engineering Division laboratory. Using a 24-hour averaging 
time these uncertainties average to about zero and the accuracy of measure- 
ments is better than 1:10". Signals at audio frequencies are also transmitted 
from the Dominion Observatory to the Radio and Electrical Engineering 
Division and from the Applied Physics Division to the Dominion Observatory, 
but these are used mainly for confirmation to lesser accuracy of measurements 
described earlier. Figures 1 and 4 show some of the results of frequency 
intercomparison. 


ABSOLUTE VALUE OF FREQUENCY 

All observations are freely exchanged amongst the three laboratories. From 
analysis of these observations coupled with astronomical observations the 
absolute frequency of the crystal oscillators is determined by one of the authors 
(S.N.K.). Figure 5 shows the results of this determination. Similar analysis 
(by M.M.T.) using shorter observational periods of 1 to 2 months gives the 
data useful in astronomy. In this work the unit of time used is a second of 
U.T.2 mentioned above. This can be later converted to a second of E.T. to 
conform with the definition of time, but in most cases this is not needed. 
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Fic. 4. Long term intercomparison of the frequency of crystal oscillators. In each case the 


ephemeris calculated from least square analysis of the observational data is listed on the top. 
The graph shows the difference between the calculated ephemeris and the actual observations. 
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Fic. 5. Absolute frequency determination by oscillator AP. The graph on the top shows 
the%frequency of oscillator AP as given by the ephemeris calculated from analysis of astrono- 
mical observations. The equation accompanying the graph is the calculated ephemeris. The 
graph at the bottom shows the difference between U.T.2 and oscillator AP used as a clock. 
Only the difference between actual observations and the calculated ephemeris is plotted. 


The accuracy of the absolute frequency determination is limited by the 
errors of astronomical observations and the uncertainty of the corrections to 
be applied. In most cases it is possible to achieve an accuracy of 5:10° if it is 
assumed that all the corrections needed for the determination of a second of 
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U.T.2 are known exactly. Recent work of Essen and Markowitz (Essen et al. 
1958) has shown that there are secular changes in U.T.2. For a further study 
of these effects and for the investigation of the possibility of redefinition of a 
time interval in terms of atomic phenomena a cesium beam frequency 
resonator has been developed by Kalra, Bailey, and Daams (1958). 


INTERNATIONAL INTERCOMPARISON OF STANDARDS 


Two techniques are used for the purpose of measuring the national standards 
of the United States of America and Great Britain in terms of the Canadian 
standard. A number of broadcast stations transmit time signals which are 
monitored by the national laboratories and observatories and their times of 
reception are published by the national observatories. The error in measuring 
the time of reception in terms of the local standard is usually about +5 x 10~¢ 
second. It is therefore seen that if the measurements are averaged over a 
month the frequency of the oscillator controlling the time signals can be 
determined to better than 1:10°. These measurements are carried out at the 
Division of Applied Physics and at the Dominion Observatory. Time signals 
transmitted by stations with the following call signals are used: CHU, WWV, 
NSS, and GBR. These observations show that the U.S. and British standards 
agree with the Canadian standard to usually within about 5:10°, which is the 
accuracy of absolute frequency determination quoted above. 

A high power broadcast station at Rugby, England, transmits at two 
frequencies, both controlled in frequency by the same crystal oscillator (Law 
1955). The transmitted frequencies are measured at the National Physical 
Laboratory, Teddington, England, against their cesium beam frequency 
standard, and the results are published in Electronic and Radio Engineer. 
These signals are compared in phase with the local standard using a tech- 
nique shown in the block diagram of Fig. 6. The rate of change of phase 
gives the relative frequency between the two oscillators. Since similar 
measurements are made at other laboratories, this technique offers a means 
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Fic. 6. Block diagram of the MSF 60-ke receiver and frequency comparator. 
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of intercomparison of standards. This work is done at the Division of Applied 
Physics and the Division of Radio and Electrical Engineering. The results 
show good agreement between the national standards of Great Britain and 
Canada. The relationship between the British cesium beam frequency resona- 
tor and U.T.2 has been described by Essen and Markowitz (Essen et al. 1958). 


CONCLUSIONS AND SUMMARY 

Over the past 3 years a frequency standard of very high precision has been 
installed in Canada. It is composed of equipment located in three different 
laboratories at Ottawa, which are separated by a few miles. Each laboratory 
has a set of quartz crystal oscillators of the Essen ring type as reference 
standards backed by working standards of different types. These oscillators 
show a stability of about 2:10" over both long and short periods of time. 
Frequency intercomparisons are all made to an accuracy of better than 1:10". 
Intercomparison of frequency between the three laboratories, which is done 
by sending signals over telephone lines, is made to an accuracy of better than 
1:10'° using 24-hour averaging time and to about 1:10° using averaging times 
of a few seconds. 

Absolute frequency is determined from astronomical observations of the 
Dominion Observatory. A Photographic Zenith Telescope and a Markowitz 
Moon Camera are used for this purpose. 

International intercomparisons of frequency standards are carried out by 
measurements of time signals and by phase measurement of low frequency 
carriers. Results show agreement with the United Kingdom and U.S.A. stan- 
dards to well within 1:10* in frequency. 
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A STUDY OF u-MESONS INCIDENT AT LARGE ZENITH ANGLES! 


B. G. WILson? 


ABSTRACT 


A Geiger counter telescope has been employed to study y-mesons incident 
at large zenith angles. The counters were hodoscoped in pairs enabling the 
direction and angle of incidence of the particles to be determined. The intensity 
has been analyzed as a function of arrival time but no significant variation has 
been found either in solar time or sidereal time. However, a variation in solar 
time has been noted on analysis of extensive air showers measured concurrently. 
The angular variation of intensity has been studied and absolute rates obtained; 
these rates have been compared with those of an earlier determination showing 
excellent agreement. Some evidence has been obtained of a directional anisotropy 
in the u-meson flux as measured at the apparatus. 


INTRODUCTION 


Systematic measurements of the intensity of cosmic rays as a function of 
arrival time have been made for many years. Several variations have been 
found relating to geomagnetic and meteorological phenomena but there is 
fairly general agreement that the radiation at some distance from the earth 
maintains a high degree of isotropy. This symmetry appears to exist even at 
high particle energies. 

To study the radiation at very high energy, experimenters have in general 
followed one of two courses. Investigation of extensive air showers at sea level 
throws light on their high energy primaries at the top of the atmosphere 
while the study of u-mesons of high energy will also yield information about 
energetic primaries. Measurements on yu-mesons have generally been per- 
formed underground to filter out unwanted low energy mesons. 

An alternative method of studying high energy u-mesons is to observe their 
incidence at the surface of the earth at large zenith angles, as suggested by 
Jakeman (1956). Since the atmosphere is not of uniform density the density 
of the air where the parents of horizontally travelling u-mesons are created 
is less than that where the parents of vertically travelling u-mesons are formed 
by a factor of about ten. Furthermore the decay interaction competition of 
the m-mesons increasingly favors decay as the path of the particle varies 
from the vertical approach to the earth. It would appear therefore that there 
is considerable advantage in making observations at large zenith angles when 
studying high energy u-mesons. 

Mesons incident at large zenith angles—‘horizontal particles’—have passed 
through 36,000 g/cm? of atmosphere so that the minimum initial energy is 
about 7X10'!°ev. The mean arrival energy at sea level for mesons- incident 


at 90° zenith angle has been computed to be 2.5X10" ev, while, at 88° it is 
1Manuscript received October 6, 1958. ; 
Contribution from the Division of Pure Physics, National Research Council, Ottawa, 


Canada. : 
Issued as N.R.C. No. 4980. a 
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1.210" ev, using the published curves of Jakeman. The horizontal flux is 
clearly distinguished by its high energy. The present experiment was set up 
to search for possible time variations in this component. 


EXPERIMENTAL ARRANGEMENT 
The apparatus consisted of four trays, each containing 16 Geiger counters, 
16 in. X1in., disposed as in Fig. 1. Counters in trays A and D had their 
axes horizontal while those in trays B and C were vertical. The telescope 





1 





Fic. 1. Experimental arrangement, side elevation. 


defined by the extreme trays A and D was 20 ft long with the maximum 
angle of acceptance 3.8 degrees above the horizontal. The telescope was 
aligned, at sea level, in the geographic north-south direction at Ottawa 
(lat. 45.4° N., long. 75.6° W.). 

Included in the telescope were 42 cm lead, while trays A, B, and C were 
shielded by 7.5 cm lead overhead and B and C by 10 cm lead at each side. 

Adjacent pairs of counters in all trays had their outputs paralleled and fed 
to hodoscope circuits and subsequently to neon lamps. The paralleled outputs 
from each tray also operated separate neons. The selection system utilized a 
Rossi-type fourfold coincidence circuit operated by the output pulses from 
the four trays. The resolving time was 4 microseconds. The ABCD master- 
pulse scanned the hodoscoped circuits operating neon lamps corresponding to 
discharged counters or trays. 

Extensive showers were detected by three unshielded trays of Geiger 
counters of total area 5200 sq. cm and one tray of area 3600 sq. cm shielded 
by 10cm lead. Each extensive tray also operated a neon lamp. Photographs 
of the neon array together with solar and sidereal clocks were taken with 
each ABCD coincidence. 

The experiment was operated from October 1956 to April 1957 with a total 
effective operating time of over 2500 hours; 11,000 events were recorded of 
which about a quarter were ‘horizontal particles’. 


METHOD OF ANALYSIS 


The events were divided into two main categories: extensive, in which any 
extensive shower tray was struck, and non-extensive. The extensive events 
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were subdivided by the selection of events in which the shielded extensive 
tray had been discharged. Extensive events therefore correspond to at least 
fivefold coincidences, the showers to contain not less than two penetrating 
particles, while ‘shielded extensives’ require at least three penetrating parti- 
cles. The non-extensive category comprises horizontal events whose directions 
of motion are uniquely determined, single counter pairs only being discharged 
in the extreme trays; events where the precise direction is not defined owing 
to multiple discharge of counters; and a small number of events which appear 
to be incompatible with single particle traversals of the telescope. 

Both the horizontal events and the extensive showers have been analyzed 
as a function of time and the significance of deviations from the mean deter- 
mined. The hourly rates have been plotted on Figs. 2, 3, and 4. 
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Fics. 2 and 3. Counting rates per hour for extensive and horizontal events as a function of 
solar time (Fig. 2) and sidereal time (Fig. 3). 








22 CANADIAN JOURNAL OF PHYSICS. VOL, 37, 1959 






R 
un 
a 





¢ SIDEREAL 


Ho 
w 
a 


? 


a oO 
2 Tet ttae tthe abt 


SOLAR 


COUNTS PER 


| 
| 
es 4 = — 


ean aesatsetea eee 
aa OSs~—CeSCSCta 5 is 2i 24 
HOURS 





Fic. 4. Counting rates per hour for shielded extensive events in scolar and sidereal times. 


Horizontal particles formed 95% of the non-extensive events and of these 
three quarters had a unique direction defined by the single counters struck 
in trays A and D. These events were analyzed as a function of their angle of 
incidence and the angular dependence of intensity obtained. After corrections 
were applied the absolute rates were obtained for each angle measured. 

The numbers of horizontal particles from the north and south have been 
compared. Finally, the production probability of secondaries has been derived. 


TIME VARIATIONS 

A total of 2296 horizontal events and 8773 extensive showers (including 
3881 shielded extensives) have been analyzed in solar and sidereal times and 
the hourly rates plotted on Figs. 2, 3, and 4. Figure 2 shows the results as a 
function of solar time for both horizontals and extensives. The statistical 
errors are also given. To determine whether the distribution of points from 
the mean is normal, the rates for 3-hour intervals were compiled and a 
x’-test was applied to the eight values to test a null hypothesis. The horizontals 
are consistent with zero variation but the extensives show a significant diver- 
gence with p = 0.001 for the null hypothesis. A sine curve was then fitted 
to the data using the least-squares method for best fit and a curve of 12-hour 
period and amplitude of 0.17 counts per hour with maxima at 0400 and 
1600 E.S.T. was obtained. This corresponds to a variation of 5.4%. 

Figure 3 shows the hourly rates plotted as a function of sidereal time. Here 
the results are not inconsistent with the null hypothesis for both cases. 

Figure 4 shows the hourly rates for the shielded extensives in solar and 
sidereal time. No established variation is apparent in sidereal time; a x?-test 
on the solar rates gives p = 0.08 for the null hypothesis, a poor correlation. A 
sine curve fitted to the data has a 12-hour period and amplitude of 0.085 
counts per hour with maxima at 0310 and 1510 E.S.T. This is an excellent fit 
with the 3-hourly averages (p = 0.75) and corresponds to a variation of 
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5.5%. Even though the shielded results are contained in the extensives, the 
apparent magnitudes of the variations and their phase angles are remarkably 
similar. 

Apart from the semidiurnal variation the results are in agreement with 
most previous work on extensive showers and u-mesons. However, the semi- 
diurnal variation noted is of interest for two reasons. Firstly, the positions of 
the maxima anc minima of the best sine wave fitted to the results correspond 
with those of the pressure wave for Ottawa. Both waves have a 12-hour 
period; the two variations, however, are 180° out of phase. If we derive a 
barometric coefficient by relating the amplitudes of the two variations, the 
results yield a value of —100% per cm Hg, much greater than the usual 
value for short-term pressure variations. 

Secondly, a variation of 12-hour period in solar time has been detected in 
high energy nucleon-induced local penetrating showers at sea level by McCusker 
et al. (1955, 1956), in high density extensive showers both in Dublin and 
Jamaica by McCusker and Reid (1958), and in very large extensive showers 
by Cranshaw and Galbraith (1957). In each case the times of the maxima 
and minima of the variations correspond closely with those of the local pressure 
variations; those of the high density showers are in phase while the other 
variations are in antiphase with the pressure waves. In these cases, however, 
barometer coefficients derived by relating amplitudes are very much greater 
than those usually associated with pressure variations of a short-term nature. 

No satisfactory explanation of the correlation and of the high barometric 
coefficients has been suggested; however, the tie-in with the pressure wave 


appears well established. 


ANGULAR VARIATION OF INTENSITY 

The 1778 horizontal particles whose directions have been uniquely deter- 
mined by comparing the single counter pairs struck in trays A and D have 
been analyzed as a function of mean zenith angle. For each of these ‘single’ 
events the direction of the particle is known within the limits of a single 
telescope. Table I shows the distribution of events; the number of available 
telescopes is given for each angle. 

Certain corrections must be made to the data in Table I to obtain an accurate 


TABLE I 
Observed numbers of events for different mean zenith angles 








Mean zenith angle No. No. Av. no. 
(degrees) telescopes events events per telescope 
86.5 2 111 55.5 
87.0 4 214 53.5 
87.5 6 226 37.7 
88.0 8 265 33.1 
88.5 10 229 22.9 
89.0 12 318 26.5 
89.5 14 280 20.0 
90.0 16 136 8.5 
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picture of the distribution. These include allowance for obstruction of the 
earth, scattering of the mesons in the lead absorber within the telescope, 
normalization for the multiple events and for particles passing outside the 
telescope but producing secondaries which discharge a tray, simulating normal 
traversals of the telescope. 

Because of the fall off of the terrain to the south of the apparatus, particles 
from the south moving upward within 0.16° of the horizontal arrive without 
interruption by the earth. Particles arriving from the north at angles greater 
than 89.7°, on the other hand, are obstructed by the Gatineau Hills so that 
the main contribution to the 90° mean zenith angle category comes from 
the south. This obstruction also affects the contribution to the 89.5° class by 
particles coming from the north but to a much smaller extent. The fractions 
of available solid angles for the different telescopes have been calculated. 

To determine the possible effects of scattering of the mesons in the lead 
absorber their arrival energies must be determined. The energy spectrum 
has been computed from the published curve of Jakeman and the mean 
energies of particles incident at 90° and 88° obtained. These energies are 
respectively 2.510" ev and 1.210" ev. 

u-Mesons of energy about 16 Gev will have a projected angle of scatter 
of 0.5°, a deflection comparable to the resolution of the telescope. It is clear 
that the average meson energy is an order of magnitude greater than that 
required to inhibit measurable effects. It is important, however, to compute 
the proportion of mesons whose energies are less than 16 Gev and evaluate 
their contribution to reduction in resolution of the telescope. 

The form of the u-meson spectrum at large zenith angles shows that low 
energy particles are strongly attenuated as a result of decay processes. At 90° 
only 10% of the incident mesons have energy less than 16 Gev and only 4% 
below 10 Gev. The number of particles which might be scattered by as much 
as the acceptance angle of the telescope is 0.5% of the total flux. At 88° 
zenith angle about 22% of the mesons have energies below 16 Gev but only 
2% may be scattered by as much as 3.8°. 

At 86.5° the proportion of low energy mesons will be greater than at the 
larger angles but particles with energy less than 3 Gev will still be limited. It 
appears likely that less than 5% of the total incident flux will be of low 
enough energy to suffer deflections comparable with the acceptance angle of the 
telescope; it is concluded that no serious error will be introduced into the 
angular distribution pattern by neglecting the effects of particle scattering. 

Multiple discharge of counters in one of the extreme trays makes it impossible 
to define the angle of trajectory of the particle within 0.5°. Multiples are 
usually caused by knock-on electrons accompanying the mesons as they 
traverse a tray and such events must be added to the single events to enable 
an absolute rate to be determined. The probability of multiple events will 
not depend strongly on the angle of incidence of the particle; particles which 
pass through the end counters of a tray will in general contribute less to 
multiple events than those through the middle counters, but such undetected 
multiples will have been accepted as singles in any case. 
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It is not possible to evaluate precisely the contribution of these undetected 
multiple events to the singles classification; however, an examination of the 
136 singles which fall into the 90° category shows that 64 passed through the 
top two and bottom two counter pairs of each tray while 72 passed through 
the telescopes defined by the central four counter pairs of the trays. Since 
the efficiency iii recognizing a multiple event is greatest for the central counters, 
then, if the contribution of undetected multiples is large, the effect should 
be greatest for the end counters. In the absence of such evidence it is con- 
cluded that the contribution of such events is small and no correction has 
been applied to the results. The multiple events have therefore been divided 
among the different angles in the same proportions as the single events. This 
results in a 27% increase in the average numbers. 

The same argument applies to corrections for particles which might pass 
through three trays and miss the fourth but give rise to secondary particles 
which discharge the fourth tray and therefore simulate single traversals. Events 
of this type are not expected to be numerous because of the small amount 
of condensed materials close to the end trays but, if they were frequent, a 
preponderance of 90° single events through the end counters over events 
through the central counters would be expected, and this is not so. The 
absolute rates have been computed and are shown in Table II. 


TABLE II 


Total number of events, effective number of telescopes, and absolute rates of particles for 
different zenith angles 
Absolute rate: 


Mean zenith angle Total no. Effective no. particles/cm?/sec/sterad 
(degrees) events telescopes 105 
86.5 141 2 7.00+0.68 
87.0 272 4 6.75+0.46 
87.5 287 6 4.75+0.31 
88.0 337 8 4.11+40.25 
88.5 291 10 2.89+0.19 
89.0 404 12 3.34+0.19 
89.5 356 12.7 2.77+40.17 
90.0 173 7.2 2.38+0.21 





The intensity of u-mesons has been plotted as a function of zenith angle 
in Fig. 5. Also plotted (open squares) are the intensities found by Jakeman. 
It will be seen that the results of the two experiments are in good agreement 
and that the values are fairly satisfied by the trend line on the figure repre- 
senting exponential absorption. 

The angular variation of intensity of the hard component for zenith angles 


up to about 60° is usually written in the form 
(1) Ig = Ip cos"6 


where J, is the intensity at zenith angle 6, Jo is the vertical intensity, and n 
has a value close to 2.0. A relation of this type, extrapolated to 90°, will 
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Fic. 5. Mu-meson intensity as a function of zenith angle. 
Fic. 6. Variation of cosine exponent ‘m’ in angular distribution relation as a function of 


zenith angle. 


yield zero intensity. However, a finite intensity at 90° would be expected 
and has been demonstrated. If equation (1) is altered to 


(2) Ig = In cos"@+A, 


where A is the intensity at 90° and is smaller than Jo by a factor of about 
three hundred, the variation of intensity over all but the largest zenith angles 
will approximate closely to that given by (1) while expression will be given 
to the known finite intensity at 90°. Substituting in (2) values derived 
from the trend line of Fig. 5 and using Greisen’s (1942) value of 0.0082 
particles/cm?/sterad/sec for J) we obtain a series of values for n, which have 
been plotted in Fig. 6. This indicates that deviations from the cosine-square 
dependence of meson intensity become significant only at the largest zenith 
angles. 
DIRECTIONAL ANISOTROPY OF HORIZONTAL PARTICLES 

The form of the angular variation of intensity shows little obvious depend- 

ence on the arrival direction of the particles. An analysis of single events is 


presented in Table III. Mesons entering the telescopes at mean zenith angle 
90° are omitted from the data since their directions cannot be specified. 


TABLE III 


Observed number of events at different zenith angles for each direction 


Zenith angle 86.5 87.0 87.5 88.0 88.5 89.0 89.5 Total 
Particles incident from north 44 72 83 113 97 121 150* 680 
Particles incident from south 67 142 143 152 132 197 147 980 


*Value normalized to seven effective telescopes. 
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There appears to be a divergence in intensity between the two directions 
of incidence. Unfortunately the results are not directly comparable since the 
mass of lead absorber within the telescope is not symmetrically placed, being 
in the southernmost quarter of its length. This might be expected to lead to 
a greater chance of multiple events being detected caused by knock-on pro- 
cesses when the particles are incident from the north, the 16-ft distance 
from the lead to the north tray biassing detection of events from the south 
in favor of singles. It is necessary therefore to analyze the multiple events 
and apportion them to the appropriate primary direction. The 442 multiple 
events were divided in the following manner. 

In 102 cases a single counter pair was discharged in tray A and two or 
more in tray D, all of which were lower in tray D than the single counter 
discharged in tray A. These events appeared most likely to have come from 
the south. Similarly 136 events appeared most likely to have come from the 
north. The remaining 204 events consisted of cases where more than one 
counter was discharged in both trays and events where the single counter 
discharged in one tray was intermediate in vertical position between the 
several discharged in the other. These events were distributed 76 to the south 
and 128 to the north on the single assumption that the number of particles 
associated with a single penetrating meson will generally increase on passage 
iirough condensed material. 

It is felt that this criterion will hold for the majority of traversals and 
that the cases in which it does not hold will balance each other to a first 
approximation. The totals have now become more even: 


ll 


particles incident from the north 680+136+128 944+31, 


ll 


particles incident from the south 980+102+ 76 1158+34. 


The difference between south and north intensities of 214 particles amounts 
to 4.6 standard deviations. 

It is important to consider experimental biasses which might be present 
because of asymmetry of the apparatus. Owing to the position of the lead 
absorber there appear two possible sources of error. First there is the chance 
of detection of secondaries produced in the lead where the primary particle 
discharges only three trays. Owing to the concentration of the lead to the 
south this would appear to favor particles from the north. There are, however, 
two balancing effects. There is a larger flux of particles through trays A, B, 
and C than through D, C, and B owing to the greater acceptance angle; also 
particles travelling through B, C, and D from the south may produce second- 
aries in the lead above tray A which could strike that tray, simulating a 
single traversal. These effects must be small since, as has been noted above, 
there is no evidence for end effects in the extreme trays when the 136 single 
events in the 90° category are studied. Further, the number of particles which 
strike the upper and lower halves of tray D has been determined for a sample 
250 events; 50.8% of the particles pass through the lower four counter pairs 
where 50.0% would be expected if there was no effect due to secondaries. 
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Secondly, the effect of scattering on the mesons will affect particles from 
the south to the greater extent as a result of the greater distance from the 
lead to the extreme tray. As has been seen, however, the number of particles 
which experience significant scattering is quite small and they are as likely 
to be scattered upwards as downwards in any case. 

Since no evidence of experimental bias has been found it is concluded that 
there is evidence of a directional anisotropy between north and south in this 
high energy radiation. It is clear, however, that there are uncertainties 
inherent in the analysis of the events which may well be as large as the 
statistical uncertainties quoted. 

If the difference in rates between north and south is not experimental in 
origin there is evidence of a north-south effect of considerable magnitude. A 
north-south effect in cosmic ray intensity is expected owing to the modifi- 
cation of the primary flux by the earth’s magnetic field and subsequently 
because the solid earth blocks many possible particle trajectories. At the 
geomagnetic latitude of Ottawa, 56.8° N., the total allowed cone is approxi- 
mately given by the combined Stérmer plus the earth’s shadow cone. Alpher 
(1950) has published curves relating the minimum arrival energy for different 
azimuths and zenith angles at 10° latitude intervals. For 90° zenith angle 
the values of Enin for particles arriving from the north and south at Ottawa 
are 14.0 and 0.5 Gev respectively for protons, i.e. the allowed cones open up 
away from the pole. However, this considerable bias towards particles from 
the south will only be apparent at low energies, measured at sea level, and 
will not affect the primaries of the u-mesons detected in this experiment. These 
primaries must have energies in excess of 10'' ev and will have a mean energy 
in excess of 10" ev. 

The anisotropy remarked in this paper cannot therefore be related to the 
influence of the earth’s magnetic field and can only be attributed to an 
asymmetry in the primary radiation at some distance from the earth. 


PRODUCTION OF SECONDARIES 


In 4588 tray traversals during the experiment, 416 events had multiple 
discharges in one tray corresponding to a probability of 9.1+0.45%. Assuming 
that the probability of discharge in any tray is independent of occurrences 
in another tray, 38 events would be expected to include multiple discharges 
in both trays. Forty-six such events were obtained. 

The average number of secondary electrons leaving a lead plate in associa- 
tion with a u-meson has been shown to be 0.07, using cloud chambers with 
horizontal lead plates, by Hazen (1943) and Brown et al. (1949). Using an 
apparatus similar to the present arrangement but arranged vertically, Jake- 
man found a production probability of 3.5%. Taking, then, an efficiency of 
50% for a counter tray in resolving and detecting secondaries, the present 
experiment shows a probability of 18+1% for the knock-on phenomenon in 
lead in association with the horizontal flux of y-mesons. With the same 
efficiency factor of 50%, Jakeman’s results yield 21+2% using a narrower 
telescope and therefore higher mean energy. Underground measurements at 
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180 meter water equivalent by Naranan et al. (1957) show the equilibrium 
number of electrons to be 17.2+1.2% using a cloud chamber technique, in 
excellent agreement with theoretical prediction. Tiffany and Hazen (1950) at 
860 meters obtained a value of 22% also using a cloud chamber. The com- 
parison of these values confirms the high energy character of the horizontally 
travelling u-mesons. 
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DIFFUSION IN MULTICOMPONENT METALLIC SYSTEMS 


IV. A GENERAL THEOREM FOR CONSTRUCTION OF MULTICOMPONENT 
SOLUTIONS FROM SOLUTIONS OF THE BINARY DIFFUSION EQUATION! 


J. S. KirKALDY 


ABSTRACT 
If C = C(x, y, z, Dt) is a solution of the binary diffusion equation 
ac 


—— ere 


for a given boundary condition, then 


n—1 


Ci = aot >, auC(x,y,2,042), 
k=1 


where 


n—1 


uy = > D 50 jx /Qix 
j= 


is a solution of 


n—1 
aC, = >. D wre. 
ot k=1 


for the z** component of an -component system subject to the same formal 
boundary condition. As an example, the ternary diffusion field about an initial 
point source is calculated. The method of calculation for steady state and 
transient equilibrium is indicated. 


INTRODUCTION 


Although a suitable phenomenological description of diffusion in multi- 
component systems has been available for some time (Onsager 1945-46; 
Darken 1951), until recently no attempt has been made to solve the equations 
and apply the results to the description of diffusing systems. This delay is no 
doubt a consequence of the natural assumption that the differential equations 
and solutions are increasingly intractable in proportion to the multiplicity 
of the system. We intend to demonstrate that this assumption is incorrect 
and that most solutions now available for binary systems can be extended 
in a straightforward manner to m-component systems. 

In recent publications, Gosting and Fujita (1956) and Kirkaldy (1957; 
papers I, II, and III 1958) have presented solutions of the multicomponent 
diffusion equation of the parametric (x/./t dependent) form for specific 
boundary conditions and for three or more components. Here we present a 
theorem which applies independently of the geometry and boundary conditions 
and which leads to a simple rule for finding the n-component solutions from 
known binary solutions. 

‘Manuscript received June 12, 1958. 
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METHOD OF SOLUTION 


Consider the differential diffusion equation for binary systems 


ac 2 
1 <= = DV’C, 
(1) a 
where D is assumed to be constant. This has a solution of the form 
(2) C = C(x,y,2,Dt). 
The generalization of equation (1) to components is (paper I) 
, ac, n—1 ‘ 
(3) a ae > Diy Ci. 
ot k=l 
where again, all D’s are assumed to be constant. We now form a trial solution, 
n—1 
(4) Ci = aiot - a .C(x,y,2,Uxt) 
k=1 


where C is given by (2) and the u, remain to be determined. By direct substi- 
tution we can show that (4) satisfies (3) provided the indicial equations 


. n—1 


(5) uy = = Diz dj /Qx 


j=1 

are simultaneously satisfied. If the a, are independent and the u, exist, this 
set uniquely determines the latter in terms of the diffusion coefficients* and 
provides a further (n—1)*—(n—1) independent relations between the ay. 
When these are combined with 2(7—1) boundary values, the (n—1)?+(n—1) 
a-coefficients are uniquely determined in terms of the D’s and the boundary 
concentration values. In general, solutions (4) will fit the same formal boundary 
conditions as (2). 

There are some limitations to the application of this procedure. For most 
physically significant solutions, the u, (as well as D) must all be real and 
positive. For m = 3 we have previously noted (paper I) that the necessary 
and sufficient conditions that «, be real and positive are that 


(6) Dy > 0 @= 2), 
and 
(7) —}(Diy— Dy»)? < DeDa < DuDe. 


The only general statement of interest in the mathematical literature con- 
cerning equation (5) applies to a symmetric matrix [D,]. In this case the 
necessary and sufficient conditions that the u, are all real and positive are 
that |D,| and all its principal minors are greater than zero. Comparison of 
this statement with inequalities (6) and (7) for = 3 suggests that perhaps 
in general it is necessary that |D,| and all its principal minors be greater 
than zero, but this is by no means sufficient. For 1 = 3 at least, the mathe- 
matical conditions are consistent with those imposed by the second law of 
thermodynamics (paper I). 


*A special case in which the u; do not exist has been considered in paper I, p. 904. The ux 
are most easily obtained by casting (5) in the matrix form to give the polynomial equation 
D—ul| = 0. The —1 solutions of this are 1), 2, .. . Un-t. 
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A difficulty with solutions of type (4) is that many apparently physically 
plausible values of the D’s and the boundary C’s lead to negative C;’s, a result 
which is quite inadmissible. If the phenomenological scheme is to be preserved, 
it must be postulated that the Dy, 7 # k, vanish with C, in a suitable fashion. 
Such a variation can be qualitatively justified on a simple statistical model 
of solid state diffusion. It has been observed for diffusion in a ternary electro- 
lytic solution by O’Donnell and Gosting (1957). The point requires further 
theoretical and experimental investigation. 


EXAMPLE OF A NON-PARAMETRIC SOLUTION 


In previous work, solutions of the multicomponent diffusion equations have 
been limited to boundary conditions which lead to parametric (x/W/t depen- 
dent) penetration curves. As an example of the application of our theorem to a 
more complex problem we consider the diffusion field about a point source 
consisting of a mixture of two solute atoms deposited at the 7 origin in a pure 
semi-infinite matrix. For the binary case, in which the number of atoms 
deposited is S, the solution of the diffusion equation in units of atoms per 
unit volume is (Jost 1952) 


oS r° 
8) | C= Gap” exp | ahah 


According to our theorem, for the ternary case 


2 2 
2 | a. ae els ae oe 
aot amex | pal exp| a 


(9) Ci = 
and 

. scot Na ah ll ee 
(10) = ant om exp| + gtam exp| 1 


If the initial point source consists of S; and S: atoms of components 1 and 2 
respectively, then 
*c 
(11) S; = 44 { Cyr’dr 
70 
and 
%00 
(12) So=4r | Cor’dr. 
ev 
Combining these relations with the indicial equations (5) we can readily com- 
plete the solutions to obtain 


ayw = Aan = 0 
ay {S\[Du—Dx+ Z)/2+Dy2S2} /8x°? D 
a2 { DS; — S2[Diu — Dig — D)/2} /82*?D 


ayn = S; ‘Sm? !? — ayy; aa = S2/82*? —an 


V { (Du — Dx)? +4D 2D} 
(DutDn+D)/2 


Il 


Il 


D 
uy; 


ig = (Di, +Dx2—D) zs 
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STEADY STATE SOLUTIONS 
Solutions of the equations 


ac, 


n—1 
(13) D> Duav’c, = —! =0 
k=1 ot 


lead from the known solutions of the binary equation 

(14) v7c =0 

in a somewhat trivial manner since, provided |Dy| ¥ 0, (13) implies 
(15) V7C, = 0. 


Accordingly, the diffusion field for a given component depends only on 
the geometry and the individual boundary conditions and is formally identical 
with that for a binary system. The multicomponent interactions appear only 
in the calculation of the diffusion currents, 


n—1 
(16) Fa = — z= DuVC,. 
k=1 


It is tacitly assumed in the above that all boundary values are fixed. It 
may, however, be that the first m of a total of n—1 components are closed 
within the system and are therefore free to seek distributions in accord with 
the constrained distributions. These free distributions finally satisfy the 
equations J; = 0, (¢ = 1, 2,...m)(de Groot 1952) while the others continue 
to satisfy (15) (k = m+1, m+2,...n—1). 


DIFFUSION THROUGH STATES OF TRANSIENT EQUILIBRIUM 


A situation may often arise in solid state multicomponent diffusion where 
some of the components have such low diffusion rates that their distributions 
can be assumed to be quasi-stationary. In systems which are finite and closed 
the diffusion fields of the rapidly diffusing components will tend to come to 
transient equilibrium with the quasi-stationary fields through the interaction 
of the cross terms in the diffusion coefficient matrix. 

For simplicity, consider a ternary system and let component 1 be rapidly 
diffusing and component 2 be slowly diffusing. We specify that the initial 
distribution of component 2 contains no infinite gradients since this would 
allow it to keep pace with component 1 in spite of its low diffusion coefficient 
(examples of this type of solution have been considered in papers I, II, and 
III). In the approach to transient equilibrium C, will satisfy the differential 


equation 

aC 
(17) a = V-DiVCi4V- DiVC; 
where 


(18) C: & C2?) 
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is the quasi-stationary distribution of component 2. When transient equili- 


brium is attained, 
(19) —J; = DyVQit+DypVCzr = 0. 


For constant coefficients and symmetrical initial conditions this gives on 


integration, 
D2 


(20) a= —— C.+A, 


where the constant A is determined by the initial condition on component 1. 
The method can be generalized in a straightforward manner to greater num- 


bers of elements. 
CONCLUSIONS 

1. With certain qualifications, the work indicates that it is now possible to 
describe analytically diffusion in multicomponent systems where solutions 
are available for the corresponding binary problem. 

2. The steady state multicomponent diffusion fields are independent and 
are determined entirely by the geometry and the boundary conditions. The 
diffusion currents are, however, interdependent. 

3. The possibility of solutions in transient equilibrium for m > 3 and the 
simplicity of their analytical treatment presages an important application of 
the multicomponent formalism to segregation problems in alloys. 
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JE™ERMINATION OF THE UNIT OF CAPACITANCE! 


A. F. DuNN 


ABSTRACT 


As part of the program to establish a wholly self-consistent group of basic 
electrical units for the country, the unit of capacitance maintained in the 
National Research Council has been determined in an absolute sense in terms of 
the frequency and resistance units maintained in the laboratories, with an 
accuracy comparable to that of the primary standards. The determination was 
made by connecting additional components to an existing capacitance and con- 
ductance measuring bridge to form a modified Wien bridge network. The results 
obtained indicate an accuracy of +0.0025%. 


INTRODUCTION 

The unit of capacitance may be derived from those of resistance and 
frequency, but of the methods which have been proposed, few yield the 
accuracy with which the standards of resistance and frequency are known 
and reproducible. The present paper describes the adaptation of a bridge 
circuit due to Wien (1891), successfully used by Ferguson and Bartlett (1928) 
of the Bell Telephone Laboratories, which permits a determination of capaci- 
tance with high accuracy. Although it is not essential for the measurement, 
the use of the Type 12 Capacitance and Conductance Bridge developed at 
the Bell Telephone Laboratories by Shackelton and Ferguson (1928) eases the 
problems of carrying out the measurement. 

In its simplest form, the bridge shown in Fig. 1 consists of two equal resist- 
ance ratio arms, a third arm containing a capacitance and resistance in series, 





Fic. 1. Basic Wien network. 
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and a fourth arm containing a capacitance and resistance in parallel. A balance 
is easily obtained by varying any two of the five variables, i.e. the two capaci- 
tances, the two resistances associated with them, and the frequency. 

The equations for balance given below are such that if the ratio of the 
capacitances and the value of the frequency are known, the magnitudes of the 
capacitances may be determined from the knowledge of only one of the 
resistances, say that in the series arm. Since the ratio of two capacitances 
may be obtained with a high degree of precision by supplementary measure- 
ments, it therefore becomes possible to use this network without a knowledge 
of the parallel resistance, the measurement of which presents certain diff- 
culties. 


THEORY OF THE CIRCUIT 


If in Fig. 1 the ratio arms are equal in resistance and time constant (or 
phase angle), then the equation of balance may be written 


(1) ee et oe 

so 

(2) | (z+) (GitjoCs) = 1, 
and by separating the real and imaginary components 
(3) Gs, = wC3C;R;, 

(4) Cs—C3+C3RsG, = 0, 


from which 
(5) wC?R; = C3/C5— a 


If the capacitances C; and C; are measured on the same capacitance bridge 
which has a wholly self-consistent calibration, then the measured values C;’ 
and C;,’ will differ from the absolute values C; and C; by a constant ratio 
K = 1+a, or 


(6) C; = KC,’ = (1+a)C,’ 
Cs = KC,’ = (1+a)C,’. 
The balance equation of the network may now be written 
ve 4 Ms 4 
(7) C= 1 | Geta) _ | = 1G 4 
wks Cs (1+a) wks Cs 
or 
| eee | 
(8) l+a = oe Ris 


Thus, the absolute value of C3, or the value of the ratio factor K = 1+a, 
may be determined from measured values of capacitance, resistance, and 
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frequency. Because many values of C3 were used, the end result of these 
measurements is a knowledge of the ratio factor K, or briefly the correction 
factor a. With this information, any subsequent measurements may be referred 
to an absolute scale of capacitance directly. 

In the foregoing, C; and Cs; were assumed to be pure capacitances and R; 
and R, pure resistances. Of course neither are obtainable in practice; the 
former will always have some slight conductance and the latter some slight 
reactance. If ca>acitors having small losses are used, the conductance of Cs 
may be considered as a resistance Rp, in parallel with Ry and that of C3; asa 
resistance Rg; in series with R3. The reactances of the resistors Ri, Re, R3, Rs 
may be included as lumped constants associated with the resistors and 
considered as a part of the over-all bridge network. 

One feature of the measurements made here is the fact that the primary 
bridge balance and all secondary balances are made with the Type 12 Capaci- 
tance and Conductance Bridge as originally designed and built. Because of 
this, it is necessary to consider the entire bridge network when analyzing any 
baiance. The schematic arrangement of the whole bridge when connected for 
either of the possible Wien network measurements is shown in Fig. 2. Both 
of the possible network arrangements shown in Fig. 2 were used, and a final 
measurement of K or a comes from the mean of the results for the NORMAL 
and the REVERSED bridges, which results in reducing the corrections due to 
the inequality of time constant of the ratio resistors from a first-order correction 
to a negligible second-order term. 

The networks shown in Fig. 2 were analyzed, and detailed balance equations 
obtained.* This leads to a value of K, or 1+a, corresponding to equation (8) 
given by 
(9) (lta)? = K* = {De Pode Or Dad 1 14 24 ae J 

Wo (D,—Dyo)s (R;*) (Dy—D4o)sRs 
for the NORMAL bridge. For the REVERSED bridge the factor in the square 
bracket becomes 
Te—T} re 
. si a acid 
[1+ ona 
In each of these equations, (D;—D4o)s and (D;—D,o)s are the settings of the 
mica capacitor dials on the Type 12 bridge when capacitors Cs and Cs are 
measured, wo is the nominal frequency of operation, R;* is the nominal value 
of resistor R3, (r2—7) is the difference in time constant of the ratio resistors, 
and Ly and Lz refer to the summation of the correction terms for the NORMAL 
or REVERSED bridges. It will be seen that the algebraic sign applying to re— 1, 
is reversed for the two bridges, and hence the term is eliminated by taking 


the mean of the results. 
% 


APPARATUS 
The basic bridge used in this determination is a completely shielded 


*Division of Applied Physics Report APEM-386. 
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Schematic diagram of Wien network on Type 12 bridge. 2a. Normal. 26. Reversed. 
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3. Wien network components connected to Type 12 bridge. 3a. Normal. 3b. Reversed. 
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equal-ratio bridge designed by the Bell Telephone Laboratories and known as 
the Type 12 Capacitance and Conductance Bridge, Style W-10105. A similar 
type of differential element bridge may be found in Harris (1952). The addi- 
tional components needed to complete the Wien network were connected to 
the bridge through the normal terminals. The schematic diagram is shown 
in Fig. 2, and photos of the complete bridge connected for either the NORMAL 
or REVERSED Wien network are shown in Fig. 3. 

In order to ensure that the additional capacitances Cep, Cceg, and C3p remain 
constant during the measurements, special care must be taken with the 
shielding of the external components. Commercially available rigid coaxial 
line elements were used wherever possible, as seen in Fig. 3, and the double 
shielded lines connecting to R; were made from special coaxial cable with 
good insulation between the shields. The resistors used for R3 were oil im- 
mersed in a metal case with a separate screen terminal, and were mounted on 
insulating pillars in a relatively large metal box which may be seen closest 
to the main bridge in Fig. 3. The inner shield of the resistor was connected to 
the junction between R; and C3, as indicated in Fig. 2, in order to specify 
capacitance C3p and keep it constant. 

The effect of this shielding has been to eliminate any detectable effect due 
to the proximity of other bodies. The stability of the shielding arrangement 
was such that the whole of the exterior network could be removed and then 
reassembled without changing any of the balances by more than one or two 
micromicrofarads; the effect on the determination of K was not more than 
one or two parts in a million. In this way it was possible to make the required 
balances with assurance. 

With any one of four resistors available for R;, there was a wide range of 
values of C;, Cs, and Ry which would satisfy the main balance equation. 
Commercially available three-dial decade capacitance boxes and a six-dial 
decade resistance box were used for these components. The capacitance or 
time constant of each of these units with their connecting lines was measured 
at the time of the determination and thus the units need to remain stable 
for a period of 30 minutes only. 

In order to establish the unit of capacitance in terms of resistance it was 
necessary to measure resistor R; in terms of known standards of resistance, 
which was done at the time of each measurement by means of an auxiliary 
Wheatstone Bridge network. The d-c. value of R; was obtained with a pre- 
cision of 4 or 5 parts in a million in terms of the basic 1-ohm standards main- 
tained in the laboratory. It was obtained within 20 minutes of the main 
bridge balances and hence the temperature stability obtained from the system 
of double shielding and oil-immersion in a temperature-controlled room is 
sufficient to warrant this accuracy. 

The power to the bridge network was derived from a crystal-controlled 
oscillator supplying a frequency of 1000 c.p.s., or from a tuning fork oscillator 
supplying 1591.55 e.p.s. (w = 10‘). In either case the frequency was monitored 
by the Standard Frequency Laboratory, and the frequency values used were 
accurate to within +1 10~° for 1000 ¢.p.s. and +10 10~° for 1591.55 c.p.s. 

The detection system used for the measurements had a sensitivity of better 
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than 0.5 microvolts at either frequency, which corresponds to detectable 
changes in balance of a few parts in 10 million for most measurements. The 
input bridge voltage was 10 volts r.m.s., at which level there were no detectable 
changes in balance due to self-heating of the components. 


PROCEDURE 


It is seen that a determination of the unit of capacitance, or of the constant 
K = 1+a, consists of the primary Wien network balance and a number of 
secondary balances. Each measurement made with the bridge consists of a 
zero balance and a measure balance with the measured value being defined by 
the difference between the two. In each of the programmed measurements, the 
zero balance was made first, followed by the measure balance, then by 
zero, by measure, by zero, and so on. Thus, if four zero balances are made, 
interspersed with three measure balances, a total of six differences are obtained, 
which may be averaged to give a mean measurement. In this way, accidental 
readings or fluctuations may be eliminated, or a drift due to temperature or 
some other cause will be seen. In the program, the measurements were made 
in the following sequence: C3, Cs, Wien, Cs, C3, 74/Rs, R3. lf there is a drift 
in value between the sets of measurements, say for C3, a mean of all the 
measurements for C; will give a result which represents the value of C; at 
the time the Wien balance was made. As C; is always smaller than C;, and 
the reading sensitivity of the bridge dials approaches the limit of error desired 
for the measurement of C;, this measurement was made as close to the Wien 
balance as possible. Altogether, for either the NORMAL or REVERSED bridge, 
there were 12 measurements of C3, 12 of Cs, and 10 of the Wien balance, as 
well as 6 of 74/R, for each individual determination of a. By taking the average 
of each group, a representative mean value is obtained, and a suggestion of 
a standard deviation or a probable error appears. 

In this determination, all the measurements are made by means of the 
Type 12 Capacitance and Conductance Bridge, using the various capacitors 
built into the bridge as standards of reference. As with all measuring bridges, 
the indications of the dials are not wholly correct and the bridge must be 
calibrated. This may be done in the present bridge by a step-down method 
of autocalibration, i.e. by comparing a step on one dial with LO steps on the 
next lower dial. Since the capacitors in each decade are completely shielded 
from those in the other decades this procedure gives an accurate comparison 
between nominally equal dial settings. When extended over the five decades 
of capacity, this furnishes a precise comparison of any bridge setting with 
any other, although it gives no information as to the absolute values of any 
of the settings. The absolute value may be determined by measuring the 
value of a known standard and adjusting the whole calibration to fit. Alter- 
natively, it may be assumed that the dial settings, when corrected for their 
relative values, differ from the absolute value by some constant ratio A = l+a, 
which is what is done in this case. It will be seen that the whole of this deter- 
mination depends on the precision of the relative values obtained trom the 
autocalibration, as well as the stability of the bridge units. 

It is believed that the autocalibrations performed on this bridge are of 
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adequate precision. In addition, comparison of results for about four years 
shows no appreciable relative shift of the values obtained. As this includes 
both air decades and mica decades of capacitance, the assumption of con- 
stancy or stability is more reasonable than postulating that two such widely 
different types of capacitor should change by the same relative amounts. 
The measurements which make up this determination of the unit of capaci- 
tance divide themselves into three distinct sets, as follows: 
Set No. 1 (December, 1956)—frequency 1000 c.p.s.; total number of 
determinations, 21. 
Set No. 2 (May, 1957)—frequency 1000 c.p.s.; total number of deter- 
minations, 20. 
Set No. 3 (June, 1957)—frequency 1591.55 c.p.s.; total number of deter- 
minations, 20. 
Set No. 1 was the first group of measurements which could be considered to 
be other than exploratory, and is given equal weight with the later measure- 
ments. As a result of the experience gained from this set, the whole circuit 
was rebuilt in a more stable form. Set No. 2 was then obtained under these 
improved conditions, but is nominally identical with set No. 1, except for 
the 6-month difference in time. Set No. 3 is a repeat of set No. 2 at a different 
frequency, but at essentially the same time. In all cases, the results were 
corrected for variation of temperature, relative humidity, and barometric 
pressure. Thus, the final results represent a value of a for atmospheric con- 
ditions of 22°C, 50% relative humidity, and barometric pressure of 29.9 in. 
of mercury. 


DISCUSSION OF RESULTS 

The values obtained from the three sets of measurements are shown in 
Tables I, JI, and III. Each value of 2a (r2—71)/CsR3 (coming from eq. 9) is 
derived from eight terms which are experimentally obtained from the main 
and auxiliary bridge balances, two terms which come from the autocalibration 
of the bridge, and one term from the temperature at which the experiment 
was made and the temperature coefficient of the capacitors in the bridge. Of 
these 11 terms which make up the value of 2 in equation 9, some are more 
precise than others, and an estimate of the standard deviation of each has 
been made. The total effect on the value of 2aF(r2—71)/CsR3 is given by 


(10) C202 = oy+o2? eee +o). 


The value of a@ itself is obtained as a mean of the NORMAL and REVERSED 
measurements, and the standard deviation of the resulting value is given by 


(2ay —2an)* 
+t 
The values of ¢ obtained by this method are shown as indeterminacies in the 
values of 2a and a in the tables. 
Examination of the data in Tables I, II, and III shows what appears to 
be a correlation between the values of a and the value of resistor R; used 


2 3 
on tor 


(11) ue = 4o.° — + 


9 


= 





R; 
(2) 


500 


1000 


2000 


5000 


(2) 
500 


1000 


2000 


5000 


C; 
(uf) 


1.030 
.790 
.610 
.450 


.810 
.590 
.450 
.350 
.270 
.160 


.620 
.300 
.170 
.080 


.140 
.120 
-100 
.090 
.080 
.070 
.060 


(uf) 


.900 
.700 
.550 
.450 
. 360 


.760 
.590 
.450 
.350 
.180 


.295 
.225 
175 
cube 
.097 


.118 
.103 
.090 
.079 


.070 


DUNN: UNIT OF CAPACITANCE 


TABLE I 


Value of a at 1000 c.p.s., December 1956 


N R 
Cs 2a—(r2—71)/CsRs 2a+(r2—171)/CsR:z a 





.012 


(uf) (X10-) (X10-) (X10-*) 
.090 —223+57 —170+60 —98+32 
.110 —177+98 —217+95 —98+49 
.130 —165+94 —240+94 —101+50 
.150 —179+98 —213+495 —98+49 
.030 —131+72 —189+73 —80+39 
.040 —139+59 —233+60 —93+38 
.050 —102+54 —191+53 —73+435 
.060 —111+52 —189+51 —75+32 
.070 —109+56 —197+52 —76+35 
.080 —37+76 —198+74 —59+53 
.010 —129+115 —166+109 —74+56 
.020 —127+83 —162+81 —72+42 
.030 —161+83 —151+82 —78+42 
-040 —130+81 —251+82 —95+51 
.007 —30+151 +19+152 —3+77 
.008 +37+143 +12+141 +12+71 
.009 —95+133 —87+127 —45+65 
.010 —111+102 —59+102 —42+53 
O11 —66+101 —75+101 —35+51 
.012 —84+95 —228 +96 —78+60 
.013 —56+90 —302+491 —89+477 
TABLE II 
Values of a at 1000 c.p.s., May 1957 
N R 

Cs 2a—(re—71)/CsR3 2a+(r2—71)/CsR3 a 
(uf) (X10-8) (1078) (X10-8) 
-100 —183+96 —359+92 —135+65 
.120 —277+95 —357 +96 —158+54 
138 —244+92 —299+101 —136+51 
.150 —184+95 —384+95 —142+69 
.158 —202+102 —275+105 —119+55 
.032 —93+71 —220+75 —78+48 
.040 —95 +62 —184+62 —70+38 
.050 —100+58 —167 +57 —68+34 
.060 —95+60 — 267 +52 —90+52 
.079 —66+70 —261+58 —82+59 
.020 —89+87 —148+80 —59+45 
.025 —90+85 —151+80 —60+44 
.030 —158+80 —144+80 —75+40 
.037 —186+75 —158+78 —86439 
.039 —105+80 —150+77 —64+41 
.008 —33+150 —142+136 —44+77 
.009 —145+127 —2+143 —382477 
.O10 —135+102 —104+103 —60+52 
.O11 —118+105 —110+109 —57+54 

—116+104 —152+102 —67+53 


a-ap 


(X10- 
+0. 


+0. 
—2. 
+0. 


—12. 
—25. 


—5. 


—7 
—8. 
+8. 


—21. 
—19. 
—25. 
—42. 


+40. 
+55. 


—1. 
“ee 
+8. 
—34. 
—45. 


aap 


(X10- 
—36.{ 


—59.¢ 
—37.§ 


—43. 


—20. 


—10. 
—2. 


—0. 
—22. 


—14. 
—6. 


—22.6 
—33.6 


—11. 


—0. 
+5. 
—16. 
—13. 
—23. 
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TABLE III 


Values of a at 1591.55 c.p.s., June 1957 


N R 


R; Cs Cs 2a—(re—71)/CsR3 2a+(r2—71)/CsR3 a a—-AaEnQ 
(Q) (uf) (uf) (10-6) (X10-8) (X10-8) (X10-*) 
500 .958 040 —93 +85 —249+96 —85+60 +13.1 
747 .050 —33+87 —232+92 —66+67 +32.1 

.600 .060 —70+87 —278+96 —87+69 +11.1 

400 .080 —101+87 —224+102 —81+57 +17.1 

. 200 .100 —40+121 —103+134 —36 +66 +62.1 

1000 .990 .010 +116+4108 —184+117 —17+94 +50.6 
479 .020 —72+74 —73+100 —36+44 +31.6 

.300 .030 —72+69 —38+88 —27+40 +40.6 

200 040 +27+482 —52+85 —6+46 +61.6 

.100 .050 —4+104 —49+101 —13+53 +54.6 

2000 304 .008 +140+91 —146+92 —2+85 +50.4 
240 .010 +62+90 —135+76 —18+65 +34.4 

.150 015 —30+87 —128+69 —39+46 +13.4 

100 .020 —44+76 —165+61 —52+46 +0.4 

050 025 —17+79 — 226+50 —61+62 —8.6 

5000 060 006 —60+88 +90+95 +7+59 +50.2 
049 007 —70+79 +52 +87 —4+52 +39.2 

.040 008 —56+79 — 103457 —40+36 +3.2 

.032 009 —96+71 — 135457 —58+34 —14.8 

020 ~—.010 —251+75 —250+69 — 125436 —81.8 








in the determination. Following a general statistical method, a correlation 
coefficient between a and G; = 1/R3 was determined, indicating strong corre- 
lation for the 1000 c.p.s. results and a questionable correlation for the 
1591.6 c.p.s. results. In a general way, this effect is to be expected from a 
discussion of the effects of intermodulation distortion in the detector.* 

The data in Tables I, II, and III show no significant differences when 
treated individually, which is to be expected if @ is to be independent of 
time and frequency over the narrow ranges considered. In the final treatment 
of the data, the three sets of values were considered as one, leading to an 
equation of correlation given by 


(12) ape = [—37—.030G;]10~6 (G3; expressed in ymho) 


with a correlation coefficient of —0.96. It is now possible to obtain a standard 
deviation for the results by comparing the individual values of @ with the 
mean values predicted by equation (12). The deviations are listed in Tables I, 
JJ, and JIJ in the columns headed (a—age). A value of a, the standard 
deviation, suggests a precision for the individual values of a group making 
up a mean value, but in itself does not suggest anything about the reliability 
of the mean value although a figure that is often quoted as representing the 
accuracy of a mean value is 3axy ba/Vn. Vhe standard deviation and 
the mean for this data amount to +£31%107® and +12X10"8 


accuracy of 
Application of the criterion that half of the residuals will be 


respectively 


*Division of Applied Physics Report APE M-386 
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less in absolute magnitude than the probable error indicates a value of 
+20X10-* for the probable error of a single observation, compared with 
+21X10-* derived from the standard deviation. This agreement indicates 
an approximate Gaussian distribution of the data obtained here. From these 
considerations, a final value of a for operating frequencies in the neighborhood 
of 1000 c.p.s. may be given as 


(13) a = [—37+15—.030G,]10-*. 


The normal use of the Type 12 Capacitance Bridge is not frequency sensi- 
tive, and the effects of intermodulation distortion are eliminated. Thus the 
appropriate value of a to be used in normal measurements is that obtained 


for small or zero Gs, i.e. 


(14) a = (—37+15)10- 


for our laboratory atmospheric conditions, i.e. 22° C, 50% R.H., and 29.9 in. 


barometric pressure. 
The figure of +15X10-° is an estimate of the magnitude of the random 


errors of the experiment. To this must be added the estimates of the various 
systematic errors which are also present, such as 


+5X10-* due to scaling the unit of resistance from 1 ohm to 10,000 ohms. 
+5X10-* representing the indeterminacy in the ratio of a-c. resistance to 
d-c. resistance. 
+1X10-® due to source frequency of 1000 c.p.s. 
or 
+10X10-® due to source frequency of 1591.6 c.p.s. 


Because the 1591.6 c.p.s. results make up only one-third of the total data, a 
figure of +25X10~* has been taken as a measure of the total error involved 
in the experiment. 

In this way, any normal capacitance measurement made with the Type 12 
Capacitance and Conductance Bridge used here may be referred to an absolute 


scale of capacitance by the relation 


Cane = ACoatats (1 +a), 
=S ACatars| 1 _ (37 +25) 10 ¢}. 


An indication of the validity of these capacitance measurements may be 
obtained by comparing results on two high quality quartz air capacitors of 
value 5000 uwuf and 10,000 wuf which had been measured at the National 
Physical Laboratory, Teddington, before delivery to Canada. The compara- 
tive measurements are shown in Table IV. 

It is to be noted that the unit of capacitance determined here is based on 
the absolute unit of resistance maintained in this laboratory by a group of 
l-ohm standard resistors, which in turn are periodically compared with 
representative units of other national standardizing laboratories through the 
services of the International Bureau of Weights and Measures in Sévres, 
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TABLE IV 


Absolute values of quartz air capacitors 


Date 5000 puf 10,000 put 


1/57 4999 .87+.02 9998. 73+ .25 

4/57 4999 .8,+ .02 Not measured 

7/57 4999 .97+ .02 9998. s3+.25 
12/57 5000.0,+ .02 9998 .8:+.05 


Note: All values reduced to atmospheric conditions of 22°C, 50% R.H., 29.9in. 
pressure, The indeterminacies listed above are the precision of measurement and 
computation, The accuracy of the measurements is +0,0025%, or +0.13 yyuf for 
the 5000 wuf unit and +0.25 pyuf for the 10,000 uuf unit. 

In September, 1956, the same capacitors were measured at the National 
Physical Laboratories, Teddington, England, before delivery to Canada. The 
values obtained at that time were 5000.04 and 9999.39 uuf. 


France. In addition, a direct comparison between our laboratory standards 
and an absolute ohm derived from our primary inductor by Romanowski 
and Olson (1957) substantiates the validity of this basis for resistance 
measurements. 


CONCLUSION 

An absolute unit of capacitance has been derived from the units of resistance 
and frequency maintained in the laboratory, by means of a modified Wien 
bridge circuit and the use of a Type 12 Capacitance and Conductance Bridge. 
It has been found that the unit of capacitance maintained by the standard 
capacitors forming part of the bridge differs from the absolute unit by 37 
parts in a million. The accuracy assigned to this determination is +0.0025%. 
This means that the measuring bridge has been constructed and adjusted 
with high precision, and for most work of a routine nature it is not necessary 
to make any correction. The unit of capacitance derived here will subsequently 
be maintained in this laboratory by a group of high quality capacitance 
standards as well as the standards in the bridge, and the bridge unit may 
be checked at any time. Periodic redeterminations of the unit of capacitance 
in an absolute sense by the method reported here will ensure that the unit 
remains consistent. 
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NOTE ON THE MOTION OF ELECTRONS AND HOLES IN 
PERTURBED LATTICES! 


R. R. HAERING 


ABSTRACT 


A new set of orthogonal, localized functions is discussed. The functions bear 
the same relation to the functions of Luttinger and Kohn as the Wannier functions 
bear to the Bloch functions. The new functions are used to give an alternative 
derivation of the effective mass equation. 


1, INTRODUCTION 


The problem of the motion of charge carriers in a perturbed periodic field 
has received a good deal of attention in recent years. For slowly varying 
potentials, Slater (1949) derived the so-called effective mass equation by 
expanding the eigenfunctions of the perturbed Hamiltonian in terms of a set 
of localized, orthogonal functions first introduced by Wannier (1937). An 
equivalent derivation can be given by using an expansion in terms of Bloch 
functions (Kohn 1958). Both of these treatments suffer the disadvantage that 
they apply only to simple band structures (i.e. nondegenerate bands with 
extremum at k = 0), and that they are not readily generalized to more com- 
plicated cases. This difficulty was overcome by Luttinger and Kohn (1955) 
through the introduction of a new set of functions. Their derivation for the 
simple band case is somewhat more involved than previous treatments; how- 
ever, it has the great advantage of being readily generalized to more compli- 
cated band structures. 

In the following sections, we shall define a new set of localized, orthogonal 
functions, and use these to give an alternative derivation of the effective 
mass equation. The new functions may be defined in terms of those of Luttinger 
and Kohn and bear the same relation to the latter as the Wannier functions 
bear to the Bloch functions. However, contrary to the derivation in terms of 
Wannier functions, the present treatment is readily generalized to more 
complicated band structures. 


2. LOCALIZED ORTHOGONAL FUNCTIONS 
We begin by defining the functions of Luttinger and Kohn (1955) for box 

normalization.” Let 

(1) ¢nk = € ” Walko- 


According to Bloch’s theorem we may write 


_ _iko-r 
Vnko =< Unko 


‘Manuscript received September 29, 1958. 
Contribution from the Department of Physics, McMaster University, Hamilton, Ontario. 


2Our definition of the functions eo differs from the definition of Luttinger and Kohn by a 


factor e~ ko't, This change is desirable so that Ko = Ynko. 
0 


Can, J. Phys. Vol. 37 (1959) 
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where 4, is a function with the lattice periodicity. Hence we may also 
write equation (1) in the form 


k k- 
(la) nk = e* “Unko 


It is readily checked that the functions ¢*¢ defined by (1) and (la) form a 
complete orthonormal set if the Bloch functions form such a set. Hence if 


(2) (Wnts Vn'k’) = Snn’ Ok’ 
then 
(3) (che, one) = Onn’ Okk’- 


Consider now the functions defined by 
k —tk-R;j k 
(4) ba'(t,Ry) =r 2d, 6 Paks 


v N“= 
= Oi Unko 2 eee, 


In equation (4), V denotes the number of cells in the crystal, R, denotes a 
lattice vector, and the sum over k extends over the first Brillouin zone. (In 
what follows, all k-sums will extend over this region.) 

Let us compare (4) with the Wannier functions, a,(r—R,), whose definition is 


1 -ik-R 
ae e z . 
J/N 2 Vnk 
A comparison of (4) and (5) shows that our 0,’s bear the same relation to 
the ¢* as the a,’s bear to the ¥,4. From the second line in equation (4), i 
is also clear that the b,’s like the a,’s are functions only of r—R,. 

The localization of 6,(r—R,) is apparent from (4). We write 


(5) a,(r—R,) = 


b, (r— R = me Unko A (¥ — R;) 


a 
VN 
where 


A(r—R,) - > ee e-Rp 
k 


A(r—R,) is a 6-like function with a spread of about one lattice spacing, peaked 
at r = R,;. Apart from a factor of (VV)-?, A(@r—R,) is in fact the Wannier 
function for free electrons. 

Our claim is now that the 8,’s form a complete orthonormal set of localized 


functions. We thus have 
(6) (bx (r—R,), Oe'(r—R,’)) = Sand jy. 


Let us first test the completeness. If we expand an arbitrary function f(r) 
in terms of the 6,’s, we write 


(7) fir) = > &,(R) oe (r—-R,). 
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Because of (4), this is equivalent to 
(7a) fe) = a Pri X) Pk 
where 


Pn (k) = aa —)>) &, (R,e™ ne, 


Hence the b,’s form a complete set if the oe do. 
The orthonormality is also readily established. From (4) we have 


(8) (ox°(r—R,), o°(r-Ry)), 
1 a 
= wey? es Rie ao ‘f aroh Cn'k’s 
= <a e™ Ri, — tk’: of ; 
NE€ k’ 
= 5 nn’d 43", 


which establishes equation (6). Finally, we point out that (4) is readily 
inverted to give 


(9) gt = 2 e* Ripko(r_R,). 


a 


The close analogy between equations (4), (8), (9) and the corresponding 
relations for Wannier functions need hardly be stressed. 


3. DERIVATION OF THE EFFECTIVE MASS EQUATION 
As an application of our localized functions we shall derive the effective 
mass equation for a simple band with extremum at k = 0. The generalization 


of this result to more complex band structures will become obvious. 
For the case of a simple band we choose the set of localized functions 


(10) b3(r—R,) = aa nod gear 


1 —tk-R; 0 
——— 2 6 
We Pnk- 


In order not to complicate the notation, we omit the superscript zero in 


what follows. 
Consider now the Schrédinger equation in the presence of a perturbation 


U(r). 
(11) Hv = Ev 
where 

H = Hot+ U(r). 


The eigenfunctions of 7 are the Bloch functions 


(12) HwWn = €n(K) nk: 
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To solve (11), we expand W in terms of the },’s 

(13) v = 2 ow (Ry )by (r—Ry). 

Substituting (13) in (11), multiplying by b%(r—R,) on the left and integrating 


over the crystal volume, one obtains a set of equations for the coefficients 
c,(R,). 


(14) 2» (nj|Ho|n'j! ow (Ry) + Dy (nj|U|n'j’)en (Ry) = Een(Ry) 
where 

(142) (nj|Ho|n'j’) = [drb¥(r—R,)Hoby (r —Ry) 

and 

(140) (nj|U|n'j’) = fdrb¥(r—R,) U(r)by (r—Ry). 


Let us look at the matrix elements (14a) and (145) more closely. Consider 
first the matrix elements of the perturbation U(r). We assume at this stage 
that U(r) is a slowly varying function of r.* Since b,(f—R,) is strongly peaked 
atr = R,;, we may then replace U(r) by U(R;) in (148), and use equation (6) 
to obtain 

(15) (nj| U|n’j’) = U(Rj)bnn5 57. 

For later convenience it is desirable to rewrite (15) somewhat. 


(15a) (nj|U|n’j’) U(R5)bnn5 53", 


a 2. Rig ®°Rs' (nk |U|n'k’), 
N ‘ek 


ll 


where 
(nk|U|n'k’) = “ bun DOM U(R)). 
Next consider the matrix element (14a). 
(16) (njloln') = 5 ee ® (ne Holn'k’) 
where 
(nk|Ho|n'k’) = [dro Hogn'x’ 


The matrix element (nk|H |n’k’) has been evaluated by Luttinger and Kohn 
(1955). One readily finds that: 


y Le s nh AK - Pan’ 
(17) (nk|Ho|n'k’) = [4 en(0) Eta, Poe Ox, 


where 
Pan’ = Jdruto(—ih¥ )up'o. (Den = 0) 


8Not all perturbations realized in practice satisfy this condition. For a further discussion, see 
Feuer (1952). 
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It is desirable to split (mj|Ho\n’j’) into a part diagonal in the n indices and 
an off-diagonal part. We write H) = H}+H} where 


1 ba et Rig &'-Ri' (nk | In'k’), 


(18a) (nj|\Ho\n'j’) = 
N €k 
(180) (nj|Hi\n'j’) = x He Rig @ Ri (nk TI n'k’), 
ick’ 
and 


(nk|He|n ’k’) = lve Nein 


(nk|H}| n'k’) = Ea 


m 


A comparison of equations (14), (15a), (18a), and (18)) then shows that the 
equations for the various c,(R,) are still coupled via the terms (j|H}\n'j’). 
We now proceed to eliminate this coupling by a canonical transformation. 
Define new coefficients d,(R,) by 


(19) cn(R;) = Ze (njle*|n'j’ dy (R,). 


The equations satisfied by the d,(R,) are then 


(20) 2» (njle~S(H0+-Hi+U )e*|n'j’)dy (Ry) = Ed,(R,). 


Here, S is to be chosen so that the equations for the d,(R,) are decoupled 
to first order. We shall not give the details of this transformation here, since 
it is closely related to that of Luttinger and Kohn (1955). If we choose S 
to have the following matrix elements in the }, representation 


(21) (nj|S\n'j!) = = ee“ ®® (nk|S\n'’), 
-) k,k’ 
where 


(nk|S\n'k’) = — (nk|Hi\n'k’), (n’ # n) 


= 0, (n’ = n) 
then the equation for d,(R,;) becomes 
(22) 7 o eR-Bi. (k)d, (Ry) +U(R;)d,(R;) = Ed, (R;). 
- J 
Here 
1a €, (0) 


é,(K) = én(0) +5 ae? k° 


is the expansion of the Bloch energy ¢,(%) about k = 0. Equation (22) ignores 
the (reduced) coupling with equations for d,(R,), (m # n). This coupling 
can be shown to be negligible if the potential is slowly varying. 








~) 


52 CANADIAN JOURNAL OF PHYSICS. VOL. 37, 1959 
Putting Rs = R,;—R, in (22) «nd using the displacement operator equation 


f(R,—Rs) = & 87 f(R,) 
we readily find 
1 +49 ,): 
23) De PRSe (k)dy(Ry)+U(Ry)dy(R,) = Eay(R,), 
which becomes, on performing the sum over Rg 
(24) lén(—7V) + U(R,)]d,(Ry) = Ed,(Ry). 


If R; is treated as a continuous variable, equation (24) is the equation 
which determines the eigenvalues E of the perturbed problem. 


4, DISCUSSION 
Equation (24) is identical with the equation derived on the basis of Wannier 
functions (Slater 1949). It would however be a trivial matter to generalize 
(24) to more complex band structures using our 6,’s, while a similar generali- 
zation is not possible in a Wannier function treatment. For instance, in the 
case of multiple minima we would simply replace the 6,’s defined by (10) by 


] si ila 
bn’ (r—Ry) = VN Unto 2 freee. 


A completely analogous derivation then leads to equation III.6 of Luttinger 
and Kohn (1955). A similar generalization is possible when the bands are 


degenerate. 
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SOME PROPERTIES 
OF THE 2.23-MEV EXCITED STATE OF P#! 


A. E. LITHERLAND, E. B. PAuL,? G. A. BARTHOLOMEW, AND H. E. Gove 


ABSTRACT 


The 2.23-Mev excited state of P*! has been studied by means of the capture 
gamma rays from the 1.70-Mev resonance in the reaction Si®°(py)P#. The 
angular correlation of the ground state gamma ray established that the resonance 
had total angular momentum 3/2, and triple correlation measurements of the 
cascading gamma rays from the compound state showed that the angular 
momentum of the 2.23-Mev state was 5/2. Coincidence measurements showed 
that the cascade gamma rays from the 2.23-Mev state to the first excited state 
at 1.27-Mev were < 3% of the transitions to the ground state. 


INTRODUCTION 

The excited states of P*! have been studied recently and the work has been 
summarized in a review article by Endt and Braams (1957). Some information, 
regarding the spins and parities of the low-lying states, has been obtained 
from the study of the correlations of the gamma rays from the reaction 
Si*°(py) P*!. The zero spin of the nucleus Si* provides an important simplifica- 
tion in the interpretation of these correlations because the only unknown 
properties of the compound state are the spin and parity. A nonzero spin 
of the bombarded nucleus introduces other parameters such as the channel 
spin mixture and the orbital angular momentum mixture which increase the 
ambiguity of a single correlation measurement (Paul et al. 1955). 

This paper is a more detailed account of the work on Si*°(py)P*, reported 
briefly by Paul et al. (1956), and also includes the results of some later measure- 
ments taken with improved apparatus. The work produced the first evidence 
that the second excited state of P*! has a spin 5/2. This conclusion has been 
recently confirmed by Broude (1958). 

APPARATUS 

The experimental equipment used for the first set of measurements, reported 
briefly by Paul et al. (1956), has been described in detail by Litherland et al. 
(1956). For the second set of measurements a conventional fast-slow coinci- 
dence arrangement (Bell et al. 1952) permitted resolving times of approxi- 
mately 50 millimicroseconds to be used. A 100-channel pulse height analyzer 
of the Hutchinson-Scarrott (1952) type* was also used in addition to the 
30-channel analyzer. 

The proton beam from the Chalk River electrostatic generator was used 
to bombard a thin enriched target of Si*®* (approximately 10 kev thick at 
1.70 Mev) deposited onto a 0.02-in. tantalum backing. 


1Manuscript received October 14, 1958. ea . 
Contribution from the Physics Division, Atomic Energy of Canada Limited, Chalk River, 


Ontario. 
Issued as A.E.C.L. No. 725. 
2Now at Physics Division, Atomic Energy Research Establishment, Harwell, England. 


*Obtained from Atomic Energy Research Establishment, Harwell, England. 
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EXPERIMENTAL PROCEDURE 


The yield of gamma rays, of energies between 3.5 Mev and 10.5 Mev, from 
the reaction Si*°(py)P*! was measured with a 5-in. diameter by 4-in. long 
Nal(TI) crystal placed 2 inches from the target at 90° to the beam. The 
yield curve, shown in Fig. 1, was studied for proton energies lying between 
1.64 Mev and 2.34 Mev. At each resonance a capture gamma ray spectrum 
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Fic. 1. The yield of gamma rays from the reaction Si*°(py)P*, with energies between 
3.5 Mev and 10.5 Mev, is shown as a function of proton bombarding energy. The arrow 
indicates the resonance studied in detail. 


was obtained. The 1.70-Mev resonance was chosen for further study because 
of the prominent ground state and second excited state gamma ray transi- 
tions. The proton energy of the resonance was determined to be 1.697+.005 
Mev by comparison with the Li’(p, ) threshold. The gamma ray pulse 
spectrum obtained at the 1.70-Mev resonance is shown in Fig. 2(a). The 
yield curve was then studied in more detail to check that only one resonance 
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Fic. 2. (a) A portion of the gamma ray pulse spectrum, taken at 90° to the beam, from 
the Si3°(py)P*! reaction at the 1.70-Mev resonance. 

(b) The low energy region of the gamma ray pulse spectrum taken in coincidence with a 
gate including the 6.71-Mev gamma rays. This figure illustrates that the predominant mode 
of decay of the 2.23-Mev state of P*! is by a gamma ray transition to the ground state. 
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was involved. Yield curve measurements at 90° to the beam of the ground 
state gamma ray showed that the background contribution at the 1.70-Mev 
resonance from other resonances was <5%. Prior to the study of the coinci- 
dence or triple correlation of the 2.23-Mev* and 6.71-Mev gamma rays, a 
coincidence yield curve of these gamma rays was taken. This showed that, 
with both 5-in. by 4-in. NaI(T1) crystals at 90° to the beam, there was <1% 
contribution to the resonance at 1.70 Mev from other resonances. It should 
be emphasized that only one target was available for the Si*°(py) P*! measure- 
ments so that it was not possible to determine whether the resonance was a 
close doublet nor was it possible to determine the natural width of the 
resonance. 

The decay scheme of the 2.23-Mev state in P*! was determined by observing 
the gamma rays in coincidence with the high energy capture gamma rays. The 
spectrum in Fig. 2(b) shows a typical pulse spectrum obtained. It is clear 
that there is no evidence for the existence of a 0.96-Mev cascade gamma ray 
transition to the 1.27-Mev excited state of P*. An upper limit of 3% was 
placed on the branching ratio. The 1.27-Mev gamma ray observed in the 
coincidence spectrum is due partly to the presence of a weak primary to the 
first excited state of P*! and to the complex cascading of higher states in P*!. 

The angular correlation of the ground state gamma ray was measured using 
a 5-in. by 4-in. Nal(T1) crystal situated with its front face 8.2 inches from the 
target spot. Care was taken to prevent the weak primary transition to the 
1.27-Mev state from influencing the measured correlation. Figure 3(a) shows 
the experimental angular correlation together with the theoretical correlations 
(Sharp et a/. 1953) for various spin assignments of the 1.70-Mev resonance. 
The least squares analysis of the angular correlation is given in Table I. The 
spin and parity of the ground state of P*! were assumed to be 1/2 and even 


TABLE I 


Angular distribution analysis 














ay S:D: a2 S.B. a3 S.D. a4 Seer: 
+0.05 0.04 —0.830 0.05 —0.03 0.06 0.003 0.05 


Note: The angular distribution of the ground state gamma ray at the 1.70-Mev resonance in Si*°(py) P®! is given 
in the table. The angular distribution was analyzed by the method of least squares into Legendre polynomials. 
The coefficients are given together with their standard deviations (S.D.). The coefficient ao is arbitrarily made 
equal to unity. 

To correct for the finite solid angle of the Nal crystal the absolute value of the az coefficient should be increased 
by approximately five per cent (Gove and Rutledge 1958). This corrected value is shown in Fig. 3(0). 


for these comparisons (Endt and Braams 1957). If the assumption is made 
that magnetic quadrupole radiation cannot compete to any appreciable 
extent with that of the electric dipole,f then it follows that the spin and 


*The measured energy of the gamma ray from the second excited state was 2.23+.03 Mev. 
All energy values quoted subsequently are from Endt and Braams (1957). 

tThe ratios 'g2/IPy1 and Ty2/l'g: for the ground state transition can be evaluated on the 
extreme single particle model (Wilkinson 1958). They are approximately 10~* and 107° respec- 
tively, which give some weight to the assumption made above. 
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Fic. 3. (a2) The angular correlation of the ground state gamma ray from the 1.70-Mev 
resonance is shown together with the theoretical curves (not corrected for geometrical effects) 
for various assignments of spin and parity to the resonance. The experimental curve requires 
a dipole-quadrupole mixture in the radiation favoring 3/2+ for the resonance. This point 
is discussed in the text. 

(b) The theoretical value of the P2(cos 6) coefficient of the angular correlation is shown 
as a function of the magnetic dipole — electric quadrupole amplitude mixture of the ground 
state transition. The experimental value of the coefficient has been corrected for the effects 
of finite geometry. 


parity of the 1.70-Mev resonance are 3/2 and even, since the theoretical 
curve for a pure dipole transition does not fit the experimental correlation. 
The measured value of the amplitude of the quadrupole radiation divided by 
the amplitude of the dipole radiation, x, can be deduced from Fig. 3(b). The 
mixing ratio x lies between —1.0 and —0.3. The large uncertainty of the 
determination arises from the insensitivity of the correlation coefficient to 
the mixing ratio in the region of «x = —0.6 and the uncertainty of the geomet- 
rical correction to the correlation coefficient. There is also uncertainty of the 
ground state angular correlation coefficient due to the <5% contributions 
from other resonances at the 1.70-Mev resonance. 

The uncertainty of the geometrical correction is due to the inclusion of 
only the high-energy portion of the pulse spectrum in the measured correla- 
tion. The calculated correction is valid only if all pulses corresponding to the 
gamma ray under study are measured as a function of angle. The correction 
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factor of approximately 0.945 (Gove and Rutledge 1958) is consequently an 
overestimate. In spite of these uncertainties the correlation demonstrates that 
dipole-quadrupole mixing occurs in the ground state gamma ray transition. 
A rough measurement of 0.1 electron volts was made of the quantity ,,[,/T 
at the 1.70-Mev resonance. If the assumption is made that I, is approximately 
equal to I, then the experimental value of T',, is small compared with the 
extreme single particle estimate of 20 electron volts for magnetic dipole 
radiation (Wilkinson) 1958. It is, however, roughly comparable with the 
extreme single particle estimate, for electric quadrupole radiation, of 0.3 
electron volts. This implies that the large quadrupole admixture in the ground 
state radiation from the 1.70-Mev resonance in the reaction Si*°(py) P*! is due 
probably to the inhibition of the dipole radiation and not the enhancement 
of the quadrupole radiation. 

Having established the properties of the 1.70-Mev resonance by studying 
the ground state transition, it is possible to study the properties of the 2.23- 
Mev state by studying the angular correlation of the primary gamma ray 
transition to that state. Experience has shown, however, that the angular 
correlations of capture radiation are ambiguous even in the favorable case 
of protons bombarding a zero-spin nucleus, provided dipole-quadrupole mix- 
tures are taken into account. The case of the triple or coincidence correlation 
is quite different. In the case to be discussed, of the coincidence correlation 
of the 2.23-Mev and 6.71-Mev gamma rays from the 1.70-Mev resonance, 
there are four unknown nuclear quantities. These are the spin and parity 
of the 2.23-Mev excited state of P*! and the dipole-quadrupole mixtures of 
the cascading gamma rays. The assumption can be made that octopole 
radiation does not play a significant role.* The theoretical triple correlation 
coefficients depend upon the geometrical arrangement of the detectors em- 
ployed for the coincidence measurements (Ferguson and Rutledge 1957). 
Consequently, by employing several different geometrical arrangements, it 
is possible to obtain information about the gamma ray cascade. In the case 
to be discussed two geometrical arrangements sufficed to determine the spin 
of the 2.23-Mev state and the dipole-quadrupole mixtures of the cascading 
radiations. The parity of the 2.23-Mev state was also assigned upon the 
basis of reasonable assumptions which are discussed below. 

Two very convenient spatial arrangements of the gamma ray counters are 
shown in Fig. 4, together with the experimental correlations of the gamma 
rays. These spatial arrangements can be called, for convenience, geometry A 
and geometry B. For the correlations shown in Fig. 4, the coincidence counting 
rate as a function of angle was measured by observing (a) the 2.23-Mev 
gamma ray in the fixed, @ = 90°, counter (geometry A) and (0) the 6.71-Mev 
gamma ray in the fixed counter (geometry B). In each case the coincidence 
counts at each angle were related to each other by counting the number of 
ground state gamma rays in the fixed counter. This procedure corrected for 


*The ratio 'y3/I'g2 can be evaluated on the extreme single particle model (Wilkinson 
1958). The ratio for the 6.71-Mev gamma ray is approximately 10~§ and for the 2.23-gamma 
ray is approximately 107°. 
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Fic. 4. The experimental results of two triple correlation measurements at the 1.70-Mev 
resonance are shown. The geometrical arrangements employed together with the gamma 
rays studied in coincidence are also shown. Curves fitted by the method of least squares are 
shown as solid or dotted lines. 


possible variations of target thickness during bombardment. Proton currents 

between 5 and 10 microamperes were used during the measurements. The 

results of a least squares analysis of the data are presented in Table II(6) and 
TABLE II 


Triple correlation analysis 








a SB. 


Geometry de S.D. 
(a) A +0.21 0.08 +0.03 0.09 
. B +0.61 0.08 —0.39 0.09 
5) A +0.34 0.05 —0.04 0.06 
( B +0.43 0.03 —0.22 0.04 





Norte: Table II gives the results of the analysis of the triple correlations of the 2.23- and 6.71-Mev gamma rays 
and the proton beam. The two different geometrical arrangements A and B are illustrated in Fig. 4. The correlations 
were analyzed into even order Legendre polynomials by the method of least squares. The coefficients are given 
together with their standard deviations (S.D.). 

The coefficients given in 2(a) were obtained from data taken at six angles with the Nal crystals at 8.5 inches 
from the target. In 2(b) the coefficients were obtained from data taken at nine angles with the Nal crystals at 


6.2 inches from the target. 
To correct the coefficients in 2(a) for the effects of finite geometry the a2 and a, coefficients should be divided 
by 0.90 and 0.71 respectively. In the case of 2(b) the corresponding numbers are 0.83 and 0.56 (Gove and Rutledge 


1958). 


are also shown as solid or dotted lines in Fig. 4. The observation of a significant 
P, (cos @) coefficient in geometry B implies that the angular momentum of 
the 2.23-Mev state is >5/2. The dotted curve shows the data for geometry B 
fitted by 1 + a2P2 (cos 0). The fit is a poor one and it can be estimated that 
if 1 + a2P2 (cos 6) were in fact the actual correlation then the experimental 
correlation would be obtained in only one in a hundred tries. Since the first 
measurement of these triple correlations briefly reported by Paul et al. (1956) 
and given in Table II(a) also gave a significant P, (cos @) coeticient in the 
geometry B correlation, the assignment of an angular momentum of >5/2 to 
the 2.23-Mev state can be considered to be a strong one. 
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Figure 5 shows the theoretical values of the triple correlation coefficients 
on the assumption that the 3/2+ resonance decays by a 6.71-Mev gamma 
ray to a 5/2+ state in P* at 2.23 Mev which in turn decays to the 1/2+ 
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Fic. 5. The theoretical values of the triple correlation coefficients are shown as a function of 
x, the magnetic dipole — electric quadrupole amplitude mixture of the 6.71-Mev radiation. 
The spin and parity of the 2.23-Mev excited state are assumed to be 5/2 and even. The experi- 
mental values of the az and a, coefficients have been increased by 10% and 25% respectively 
to correct for the effect of finite geometry. After this figure had been prepared more accurate 
estimates of the corrections (Gove and Rutledge 1958) indicated that they are probably 
underestimates. The conclusions in the text are not altered by the change required in the 
figure. 


ground state. These theoretical coefficients were obtained from Ferguson and 
Rutledge (1957). Good agreement between theory and the experimental 
coefficients taken from Table I1(d) for this sequence is obtained for values 
of the mixing parameter x lying between —0.3 and —0.15.* The agreement 
between theory and experiment for the less accurate coefficients taken from 
Table II(a) is less satisfactory because the values of the mixing parameter x 
needed to fit the correlations in the two geometries are different. This dis- 
crepancy provided one of the reasons for repeating the measurements. The 
measurements, however, support the assignment of 5/2+ to the 2.23-Mev 
state in P#, 

An angular momentum assignment of 7/2 is completely eliminated by the 
data. This elimination can be illustrated most clearly by considering the 
ratios of the az coefficients obtained from the two geometrical arrangements 


*This conclusion is not sensitive to the uncertainties in the correction to be applied to the 
correlation coefficient due to the finite angular spread of the counter. The maximum values 
of this correction are given in Table II. 
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employed in the measurements. Since the same distances of the counters 
from the target were used for the geometry A correlation and for the geometry 
B correlation the ratio of the a2 coefficients for the two geometries is indepen- 
dent of the distance of the counters from the target. The values of this ratio 
obtained from Table II(a) and II(0) are 0.34 + 0.10 and 0.79 + 0.13 respec- 
tively. The value calculated for an intermediate state spin of 7/2 is —0.26. 

A useful feature of the triple correlation method can be illustrated by 
considering as an alternative to the above assignment the following possi- 
bility. The 2.23-Mev gamma ray need not be a secondary transition since 
the time sequence of the 6.71-Mev and 2.23-Mev gamma rays is not deter- 
mined by the apparatus. Let us suppose, then, that the 2.23-Mev gamma 
ray is a primary to a state of spin and parity 3/2+ at 6.71 Mev in P*. For 
this cascade the geometrical arrangements A and B, Fig. 4, are interchanged. 
The original geometry A becomes the new geometry B’ and the original 
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(b) GEOMETRY B’ 


Fic. 6. The theoretical values of the triple correlation coefficients, for the geometries 
A’ and B’, are shown as a function of x, the dipole-quadrupole amplitude mixture of the 
6.71-Mev radiation, and y, the dipole-quadrupole amplitude mixture of the 2.23-Mev radiation. 
The spin and parity of the resonance and the 2.23-Mev state are assumed to be 3/2 and even. 
= results of the calculations are shown as contour plots of the correlation coefficients a: 
and a4. 

(a) The a2 and a, coefficients for geometry A’. (b) The az coefficient for geometry B’. 
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geometry B becomes the new geometry A’. It is possible in this case to have 
a nonzero d, coefficient for the geometry A’ correlation even though all spins 
involved in the gamma ray cascade are <3/2. The theoretical correlation 
coefficients are shown in Fig. 6. It is possible to obtain values of the mixing 
parameters x and y from the experimental geometry A’ correlation coeffi- 
cients. However, the theoretical geometry B’ correlation deduced from the 
knowledge of x and y is then very badly in disagreement with experiment 
(a2 ~ 1) so that the experimental correlations do not agree with the above 
assumptions. By similar arguments it can be shown that the experimental 
correlations also do not support the possibility that the 2.23-Mev and the 
6.71-Mev cascade proceeds through a 5/2+ state at 6.71 Mev. 


CONCLUSIONS 

The 2.23-Mev state in P*' has been shown to have an angular momentum 
of 5/2 by exploiting the technique of triple correlations. The parity of this 
state is undetermined unless the assumption is made that magnetic quadrupole 
radiation cannot compete significantly with electric dipole radiation. If this 
assumption is made, then the parity of the 2.23-Mev state is even since the 
transition from the 3/2+ capturing state is experimentally a mixed dipole- 
quadrupole transition. 

In conclusion it is interesting to point out the similarities between the 
first three states of Si? and P*. The available information (Bromley e¢ al. 
1957), and the work reported in this paper, are shown in Fig. 7. The strong 
E2 crossover in each case has stimulated the application of the collective 
model to Si” (Bromley et a/. 1957) and P* (Broude, private communication). 
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Fic. 7. A comparison of the experimental properties of the low states of the nuclei Si” 
and P#!, 
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ON THE GENERAL PLANE PROBLEM OF PLASTICITY 
AND ITS GEOPHYSICAL SIGNIFICANCE! 


S1str CHANDRA Das? 


ABSTRACT 


The results obtained from the fault plane studies made in the Dominion 
Observatory show strong horizontal displacements in the focus of earthquakes. 
This indicates a marked deviation from the existing theories of geotectonics. To 
account mathematically for them, the study of the plane problem of plasticity is 
proposed using as yield condition a general functional relation between the 


stresses. 
The differential equations involved in the problem are of non-linear type. For 


solutions these are replaced by a different set of completely equivalent equations 

expressing variations along families of curves, known as characteristics, across 

peste certain derivatives may be discontinuous under suitable boundary con- 
itions. 

Physically, the elastic-plastic (or rigid-plastic) interface is taken as a character- 
istic curve which is a fundamentally unknown thing, and which must be ascer- 
tained from symmetry of the problem or by proper experimentation. Investi- 
gations in this direction are now being done in connection with the International 
Geophysical Year. Meanwhile here we discuss the various relatively modern 
methods of linearization of the stress equations. Each of them essentially deals 
with transformations of co-ordinates, depending upon the form of the yield 
condition, though in every case the generalized nature is maintained. 


INTRODUCTION 


One of the important problems facing geologists is to explain the orogenic 
belt which surrounds the Pacific basin. This belt is made up of a number of 
mountain ranges, island arcs, and ocean deeps, arranged parallel to the 
boundary of the basin, usually with associated volcanoes. That the zone is 
active is proved by the great number of large earthquakes associated with 
it. These earthquakes have foci at depths ranging from the surface to more 
than 600 km, and appear to define zones, approximately planar, striking 
parallel to the basin and dipping beneath the continents at angles of 45° 
to 60°. 

Most explanations of this orogenic belt have suggested that it is the result 
of forces acting normal to the Pacific boundary; this would explain the folding 
up of mountains and the downthrusting of deeps parallel to the boundary. 
Because the earthquake foci appear to define a plane, it was assumed that 
this was the fault plane of the earthquakes and that motion was either up 
or down this plane. Such faulting would have a strong vertical component 
which is consistent with the vertical displacement in mountains and deeps. In 
recent years seismologists (Hodgson 1957; Keylis-Borok 1956) have been 
investigating the nature of this displacement and find that the usual assump- 
tions are not valid. Most of the faults are nearly vertical, the planes do not, 
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in general, coincide with the plane of the foci, and the displacement is almost 
purely horizontal with only a minor vertical component. 

How is this apparent discrepancy to be reconciled? It may be that the forces 
do not act normal to the ocean boundary, or it may be that the properties of 
the material are different from those usually assumed. Because there are so 
many reasons for supposing that the forces are normal to the Pacific boundary, 
the second possibility appears to offer the more fruitful field for investigation. 

Most geological studies regard the earth as perfectly elastic but, since we 
are dealing with conditions at depths as great as 600 km, where the tempera- 
tures and pressures are very high, this is not satisfactory in the present 
problem. Nor can we regard the material as viscous; viscous material cannot 
fail by fracture, and earthquakes are certainly the result of fracture. We 
conclude that the material must be regarded as plastic, and the problem 
becomes one of defining the yield conditions of this plastic material. 

Much of the literature on plasticity, particularly that relating to the plastic 
failure of metals, has made use of particular yield conditions, either those 
defined by von Mises, or those defined by Tresca. These two conditions 
become the same in two-dimensional problems, and they have the advantage 
of much physical confirmation in the field of metals. However, it seems doubtful 
if they can properly be applied to the present problem, considering the very 
high temperatures and pressures involved, and the fact that these are active 
for a very long time. One might search for other specific yield conditions but 
it is difficult to justify any particular set of conditions in the absence of 
physical evidence. Moreover there is no reason to suppose that the same 
yield conditions operate at all times or at all places. It will therefore be 
informative to discuss the problem in full generality, taking the yield con- 
ditions in the form of a general functional relationship. That is the purpose 
of the present paper. 

It should be mentioned that a considerable literatuve exists in a field 
intermediate between that of the failure of metals and that of the failure of 
the earth. This is the discipline of soil mechanics in which the soil is usually 
replaced by an idealized material which behaves elastically up to some state 
of stress at which slip or yielding occurs. A certain amount of similarity 
between Prandtl’s theory of plasticity (1920) and the theory of earthy media 
has been shown by Sokolovsky (1939). In more complete plane investigations 
an extended Coulomb’s rule is used (Terzaghi 1943). Implications of assuming 
the soil to be a perfectly plastic body have been shown to be far reaching by 
Drucker and Prager (1952). The former author (Drucker 1953) has further 
extended the idea of limit analysis towards an interesting development of 
complete solutions. Thus it is necessary in such a situation to study general 
yield conditions; this supports the conclusion that general yield conditions 
are necessary in the larger problem of the failure of the earth. 

It must be understood that the investigation of the properties of the 
materials of the earth’s interior is much more complicated than the equivalent 
problem in soil mechanics because data on the properties of the materials are 
seldom known conclusively and laboratory experiments are not possible. 
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Further, the usual assumptions of incompressibility and isotropy of elasticity, 
and the usual plasticity theories, most likely cannot be taken for granted in 
our case. For example, the stress-strain relations in horizontal directions and 
in vertical directions should not be taken to be the same, and change in 
volume is naturally expected under high pressures and temperatures. 

The present analysis is limited to the plane problem. This simplifies the 
mathematics, and, as in other branches of applied mathematics, we may 
expect it to produce useful results through the proper orientation of co- 
ordinates. Indeed it has already done so; the two-dimensional study of 
approaching parallel rigid plates has shown that, for viscous fluids or for 
plastic or Bingham substances, flow can take place along directions per- 
pendicular to the applied forces (Das 1958). This helps to explain the observa- 
tions mentioned earlier that the motions in large earthquakes may be transverse 
to the forces causing them. 

The practical limitations of these studies must be emphasized. We have 
no knowledge of the yield point or the stress-strain relationship appropriate 
to any point or zone, although observations are being made in this direction 
in both the U.S.A. and the U.S.S.R. in connection with the I.G.Y. (Gzovskii 
1957; Risznichenko 1957). Moreover, there are not very many problems for 
which the boundary values can be stated in such a way as to ensure the 
existence and uniqueness of solutions. Nevertheless, the knowledge obtained 
of the pattern of stress state and of the process by which fracture develops 
will allow a better understanding of the physical situation. This in turn will 
allow a more precise evaluation of the problem. 

The development of the solutions of the general plane ‘problem is com- 
paratively modern and has not yet been included in any monogram of plasticity. 
The problem is rather difficult. It is non-linear and, moreover, the solutions 
are not usually orthogonal nor everywhere hyperbolic. The determination of 
the regions in which the solutions are not hyperbolic must be made at the 
same time that the solutions are obtained. 

Three methods of solution will be discussed in this paper, those of Neuber 
(1948), Sauer (1949), and Sokolovsky (1955). Geiringer (1953) has already 
reviewed the work of Neuber and Sauer on the linearization of the stress 
equations, but apparently was unaware of an earlier paper of Neuber (1947) 
in which he gives a somewhat fuller treatment of the anisotropic case. 

In Neuber’s method for linearization of the stress equations, (¢;—¢,) and 
t,, have been used as stress variables. Sauer takes o, and 1,, instead, and 
derives a nonlinear differential equation of second order, which he then 
transforms into a linear equation by means of Legendre transformation. Both 
methods are capable of applications in anisotropic medium as well, and equally 
simple construction of the characteristics have been obtained in both cases. 

Sokolovsky extends the usual method of solution for von Mises (and 
Tresca's) condition of yield to cover the general condition of yield expressed 
as a function of mean normal stress and shear stress intensities, without 
assuming the incompressibility of the material. Different methods of trans- 
formations have been indicated to deal with both hyperbolic and elliptic 
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nature of the equations of characteristics and solutions have been obtained 
in the form of trigonometric series. 
METHODS OF LINEARIZATION 


1. (a,—ay) and rz, as Stress Variables 
The usual equations of equilibrium in plane problem in cartesian co-ordinates 


are 
9 Feo 
(1.1) Ox oy 
Or Oa 
a a 
Ox - oy ®. 


where oz, a, are the normal components of the stress and 7,, is the shearing 


stress. 
The compatibility conditions of elasticity are replaced by the condition 
of plasticity or the law of plastic flow. This is a functional relationship between 


the stress components and is written in its most general form as 


(1.2) f(x, oy, Try) = 0. 
It is to be noted that this relationship is applicable to both isotropic and 
anisotropic cases, i.e. it is not necessarily composed from the invariants. 

All attempts to solve these equations consist of different ways of trans- 
formation to different systems of variables. 

The basic equation is obtained by Neuber in two stress variables ¢,—c, 
and 27,, in a system of oblique co-ordinates u, v. Denoting the inclination 
of u and v directions with the x-axis by a, and a, respectively (Fig. 1), equations 
(1.1) take the form 


| ee ae ra} 3 aes 2a] 
AI sina a COS ay ree i SIN Qy 3 COS ay, av = 0, 


a a Bee Ua Oty oe. 
AY sin Oo COS ay x sin ay, ay COS ay ay = (0, 


where h,du and h,dv are the line elements. 
Restrictions are now imposed upon a, and a, in such a way that the equa- 
tions (1.3) are obtained in terms of derivatives with respect to u only in the 


first equation and v only in the second by taking 


(1.3) 








; 00 Or 
sin a, —"—cos a, —" = 0, 
(1.4) ov ov 
; Try Oo, 
sin a, —=~"—cos a, =" = 0. 
” Ou ” Ou 
Then equations (1.3) become 
‘ Oo, OTzy 
sin a, ——— cos a, —— = 0, 
(1.5) Ou Ou 
; Or 0a 
sin a, —"“—cos a, —" = 0. 


Ov Ov 
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Vv 
(u = CONST.) 







(Vv = CONST.) 





From these two sets of equations we get 


2 
We. de4_( Star) =0 along v = const. 
Ou Ou ou 


Oo, Io, _ (2) : 
ov ov ov 


Therefore for both these families of curves in general we have 


(1.6) 
0 along “ = const. 


(1.7) (do,)(do,) — (dr;,)? = 0, 


which can be written in the form 
(1.8) [d(o,+0,)]? = [d(o,—¢,)]?+[d(27,,)]?. 


The yield condition expressed in terms of ¢,+¢,, ¢,—o,, 27, can now be 
interpreted as a surface called a ‘‘plasticity surface” (Fig. 2). 

Equation (1.8) clearly indicates that the characteristics of the differential 
equations of the general two-dimensional plasticity problem correspond to 
the 45° slopes of the plasticity surface. On the plasticity surface the lines 


“ = constant or v = constant constitute 45° slope lines for the angle a (i.e. 


a, OF @,). 
From equations (1.4) and (1.5) it clearly follows that 


tan a = do,/drz 


and 

(1.9) cot a = do,/drzy. 

Hence 

(1.10) cot 2a = d(a,—«,)/d(2rzy), 


i.e. 2a is the angle between the tangent to the projection of the line v = constant 
(or u = constant) on the (¢,—«,), (27,,)-plane, and the (¢,—¢,)-axis. 
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u 
Vv = CONST. 





Fic. 2, The line v = const. on the plasticity surface and in the drawing plane (after 
Neuber). 


The integration now can be done by the method of projection. Substituting 
§, n for o,—o, and 27,, it is observed that x, y-plane is projected to &, 7-plane 
by transforming the wu, v-lines of the x, y-plane into the 45° slope line of the 
plasticity surface projected on to &, n-plane. The 45° slope lines give a unique 
picture of the x, y-plane and its state of stress. 

It should be noticed that the inclination of the tangent to a line u = const. 
or v = const. is equal to half of the inclination of a tangent to a line v = const. 
or u = const. respectively, in the corresponding point of x, y-plane. 

The use of the stress variables ¢,—¢, and 2r7,, (i.e. &, 7-plane to be the 
basic plane of the plasticity surface) has the advantage that slope lines can 
be very easily constructed from the contour line (¢,+¢,) = p = const. 


2. o, and T,, as Stress Variables 

Sauer uses ¢,;, Tz, as stress variables. His method also covers the anisotropic 
case. He deduces a second order nonlinear differential equation and then 
transforms it into a linear equation by means of Legendre’s transformation. 


The flow condition is written in the form 


(2.1) Oy = f(z, Try): 
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The potential function ¢(x, y) defined by 


—— a. 


(2.2) Cz ay ’ Try = dee 


is now introduced and the quasi-linear equation of second order, 


a*6_af ao, of a6 _ 4 


2. = 
(23) ax” ax ay | atz, Oxdy 
is obtained. Applying Legendre’s transformation 
0 
uu = 5 = —Try 
0 
v= = = Or, 
(2.4) & = ux+vy— = 0, ¥—T1 X— $, 
O® o® 
s= = -— 
ou OT ry 
o® Oo® 


oe oF ao 
. Ov do,’ 
in (2.4) we get the homogeneous linear equation, 


vo af ob, Hf _dr® 





= 0, 


(2.5) 3— iT. 
OGz 00s Oxy Ofzy OGeOTcy 





where 
Ou, W , W, Ww 
dx dy” Ox dy’ 
which we assume to be the case. 

The function f(¢,, 7,,) as well as its derivatives df/d0, and df/drz, may 
be quite complicated when the yield condition is symmetric. For the charac- 
teristics in the physical plane from (2.4) we find 
; 2 Of ,2 of 
(2.6) dy —~— dx" —~-— dy dx = 0 

Oo, OF xy 


and in view of u = —r,, and v = o, the compatibility conditions 
(dex) es (=) 

try) 1 dy/ »’ 
(%.) 7 (#) 

dtz,/ 2 dy/ 


are obtained where 1, 2 denote two characteristics. 
From this, both corresponding families of characteristics of equations (2.6) 


(2.7) 


are obtained as 


0d 9 of ° 
of do, — oY do,dt,, = dt,, —do,.do, = 0, 
0c, OT zy 


(2.8) dtr 


which coincides with (1.7). 
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The characteristics are the slip lines of the present plastic stress field and 
from their properties methods have been developed by Sauer to construct 
them. 

The slip line nets like the o,,7-characteristic nets corresponding to them 
consist of two orthogonal families of curves when df/dc, = 1. They are thus 
characterized by the form 


(2.9) Or— Oy = Q(Tzy) 


where g(rz,) is an arbitrary function of 7;,. 
If this is the case we get from (2.9) 


ay (Be), = a st o/ 44). 


The o,, 7 characteristics intersect the lines 7,, = const. at a fixed angle 
and are therefore parallel and congruent. 

It may be remarked in this connection that the orthogonal family of nets 
will be particularly useful in explaining the probable fault planes obtained by 
geometrical methods from the recordings of a single earthquake made from 
different places. Taking for granted that fault movements occur along slip 
lines, which most likely is the case, they are obtained as along either of two 
perpendicular planes. Because of the uniqueness of null vectors (Hodgson 
1957) each of them denotes some point, in the plane case, where the two 
slip lines, each of different families, meet. 





3. Extension of the Classical Method 
The general condition of plasticity has been written by Sokolovsky in 


the form 
P(o,7) = 0, 
or 
(3.1) r = 7(0), 
where o and 7 are the mean normal stress and the intensity of stress re- 
spectively. 
Introducing s and ¢ defined by 
= 2(or+oy,) 
and 
= V4 (G2— oy)? +7247}, 
the conditions of plasticity for both plane strain and plane stress can be 
written in the form 


F(s, t) = 0, 


or 
(3.2) t = t(s). 


Taking F as a plastic potential and following Drucker and Prager (1952) 
we can have 
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€x ey €:  2Yzy 
(3.3) dF OF OF OF 


Oo, Oo, Oo, Ory 














€r,- ++» Yry--+ being the components of strain. Introducing e and g, defined by 
e = 2(eite2) = dete), 
g 2 (€1— 2) V {3 (Ge €y)?°-+ 79") ’ 


é, and é2 being the principal strain (e; > e2), we get 


ll 
ll 


€,—€ = g/t(or—S), 


(3.4) €y—€ = g/t(o,—S), 
é—€ = g/t(o,—S), 
Yr = g/trey 
and 
(3.5) e = g{(dF/ds)/(AF/dt)} = gh 


where / is a known function of s and ¢. 
The stress components can also be written in the form 


or = S+k cos 2¢, 
(3.6) oy, = s—tcos 2¢, 


t sin 2q, 


Try 
where ¢ is the angle between the x-axis and the direction of principal stress. 
And similarly for strain components we get 

€, = e+g cos 2¢, 

(3.7) €, = e—g cos 29, 

Yry = g sin 29. 


Substituting the values of ¢;, oy, tzy from (3.6) in equations of equilibrium 
(1.1) we get 


i Oi a 28) 
ax 1 OOS 2 ay t Sin 2¢ ay 2t (sin 2¢ ax COS 2¢ ay 0, 
(3.8) 
ds, . at dk do, . a8) 
malt PR As a gegen OM capa Oa ee hoa) 
ne 26 ax 08 2¢ eo (cos 2 ax tin 2 ay 0. 


The procedure is exactly the same as that followed in the case of classical 
plane problems with Tresca’s (and that of Mises) yield condition. The only 
change in equation (3.8) is due to the fact that ¢ is not a constant in this case. 

Further from (3.7) we can have 


. Ou ov , Ou 
Piet Ore Seen pcs i es eer 
wil 2 sin 2¢ - (h+cos 29)( 2424) 0, 


2 sin 2 ~— (h—cos 29)( 42H) 00 
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The characteristics of equations (3.8) and (3.9) are of the form 


2, 
doe V =) gins =0 





2h , 
1 dy _ sin2¢#/(1—h’) 
(3.10) dx cos 2¢—h . 
dv _ sin 26++/(1—h’) 
du cos 2¢+h , 


It is clearly seen that these equations are not confined to the hyperbolic 
case alone but may well be elliptic in type depending on whether |h| < 1 
or |h| > 1. 

In the hyperbolic case (|k! < 1) introducing y and \ defined by 

dt/ds = —h = cos 2y 
and 
3.11) 2dd = +/(1—h?)(ds/t), 
sin 2y ds/t, 
tan 2y d Int, 


ll 


I] 


the characteristies of the system of equations (3.8) and (3.9) are determined 
by the equations 
n = A(~)—¢@ = const., 
(3.12 dy/dx = tan (@—y), 
du/dv = —tan (@¢—y); 
— = A(v¥)+¢ = const., 


dy/dx = tan (@+Yy), 
du/dv = —tan (@+y). 


Il 


~~ 
_ 
w 


Jt should be noted that characteristics in the x, y plane are always orthogonal 
to the corresponding characteristics in the u, v plane as from (3.12) and (3.13) 


it is seen that for both families 


at (dy/dx)(dv/du) = —1. 
By further substitutions given by 

Z = - [x sin (+ p)—y cos (9+ y)] 
sin 2y 

q = 4 os ( - _) Si ( a 

9 = 9g ly om (O-¥)—x sin (-¥)] 

a lu cos (¢—y)+vu sin (d— p)] 
sin Zy 

§ = {fu cos (@+yp)+vsin (d+ yp)] 


~ sin zy 
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and then introducing new variables, X, Y, U, V, defined by 


/ 





Bs Bal an 
sa" (t sin 2y, 
G8 eas 
> es tp a. 
and finally putting 
if 
PQ) = sin Qy ° 


equations (3.12) and (3.13) become 
(3.14) n=XMY)-¢ 
(3.15) é=(¥)+¢ = const, dX+YVPd\ =0, dV = VPdy, 


from which the canonical forms for co-ordinates and the velocities can be 


const., dy+XPdi =0, dU = VPad,, 


I 


obtained as 


aY Py at Be 
(3.16) “ae 9 A = 0, an 2 } 
and 

ne ou P- OV P : 
(3.17) OE = 2 s On  - l 


respectively. 

Several particular solutions, namely, for the cases when (i) the family of 
characteristics (8.15) consists of straight lines £ = const., and when (ii) the 
family of characteristics (3.14), in the (x, y) plane, consists of straight lines 
n = const., and finally when (iil) both the families of characteristics (3.15) 
and (8.14) represent systems of parallel straight lines, have been discussed 
by Sokolovsky. 

Sokolovsky has also dealt with another method of transformation where 
he has extended his method (Sokolovsky 1950) for r = & (const.) to suit the 
general form of plasticity condition. Following the practice in plane gas flow 
he has obtained the results in the form of trigonometric series. The method 


is applicable for both hyperbolic and elliptic cases. 


SUMMARY 


The different solutions of the problems discussed in the paper cannot be 
compared with observation because of the generality of the assumptions, but 
they are believed to be of geophysical significance, particularly tor explaining 
deep focus earthquake mechanism and the probable fault movements there 
The geometry of the methods discussed here can be used to determine the 


flow (or fracture) pattern, if the boundary stresses as well as the particular 


yield condition are known, Conversely it will also help to understand the 
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nature of stresses and the yield condition provided the displacements are 
ascertained. It is both interesting and useful to proceed with various yield 
conditions and to find their influences on the solutions. The conclusions can 
only be decisive when comparisons with actual observations and experiences 


will be possible. 
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AN ISOTOPE EFFECT IN ANTIFERROMAGNETISM AT LOW TEMPERATURES 
J. M. DANIELS 


It is known, as a consequence of Nernst’s heat theorem, that all para- 
magnetic substances must become either ferromagnetic or antiferromagnetic 
at low enough temperatures. In most of the substances so far investigated 
experimentally, the low temperature state appears to be antiferromagnetic. 
The Néel temperatures of many hydrated salts containing paramagnetic ions 
are well below 1° K, for example KCr(SO,4)2.12H2O has a Néel temperature 
of about 0.01° K. The purpose of this note is to point out that the macro- 
scopic properties of the antiferromagnetic state can be greatly influenced by 
the isotopic composition of the paramagnetic ions, because the hyperfine 
coupling between the unpaired electrons and the nuclear spins of the para- 
magnetic ions can be much larger than the coupling between individual 
ions. As examples, let us consider Mo*+ as an impurity in K;InCl..2H,0, 
and also Sm(C2H;SO,);.9H2O. These examples are chosen because these ions 
have an isotropic g-factor but an anisotropic hyperfine coupling. 

The Mo** ions can be represented approximately by a spin-Hamiltonian 


(1) KH = gBH.S8+ AS,I1,+B(S,I,+S,1,) 
where g =~ 1.95, A =~ 0.008 cm, B =~ 0.004 cm", 


8 is the Bohr magneton, and J is the nuclear spin; in this case J = 5/2 
for the odd isotopes Mo” and Mo. (For a definition and explanation of the 
spin-Hamiltonian, see Bowers and Owen (1955), where the figures quoted are 
tabulated.) It is clear that for the even molybdenum isotopes the hyperfine 
interaction is absent (J = 0) and the antiferromagnetism is that of free 
isotropic spins; hence the anisotropy energy is expected to be small. For the 
odd isotopes, the anisotropic hyperfine coupling introduces an anisotropy 
energy, and we now derive an expression for this. 

Let us consider just one sublattice of spins all with the same crystalline 
field (s) axis. Assume that we can replace the interaction between spins by 
an effective field JZ whose magnitude is independent of direction. In a typical 
paramagnetic substance (though perhaps not for this, for which few data 
exist), JZ would be about 250 gauss. The total splitting in this field is 
g8H ~ 0.063 cm, while the hyperfine splitting amounts to 5A/2 or about 
0.02 cm~!. Under these circumstances, one can find the energy levels of (1), 
treating the hyperfine interaction as a perturbation (see, for example, Bleaney 
(1951)). It is then a straightforward though tedious matter to write down 
the free energy F = kT In Ye~*7/*7, At the absolute zero, however, F is equal 
to the lowest eigenvalue of AW; this is to second order 
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_ say SK 28 vecin's (BO4’)’ L 
gbBH 4 16 08 2 sin 8 K 2¢8H+ (5K /2) 





BS a(4 =) 1 
16 K 2g8H-+ (3K /2) 
where K? = A? cos*@+B? sin’, and @ is the angle between H and the z-axis. 
Expanding this, we have 


mr 42_ pR?2 ; 

(2) Fy-+c08%0 e a 5 +0(-45) |+-0%coe'a 

where Fp is that part of the free energy which ducs not contain @. Now @ is 
the angle between z and H while the anisotropy energy is usually determined 
in terms of an angle ¢ between z and S. The components of the magnetic 
moment vector », which is parallel to S, can be found from (2) using 
uz, = OF/0H,; u, = AF/0H,, from which it follows that the difference between 
6 and ¢ is negligible if A, B < g8H and @ and ¢ are small. According to one 
convention for defining H, the anisotropy field (Nagamiya, Yosida, and 
Kubo 1955), we have from (2) 


f5 A*—B" 
\8 B 
Thus, in the case of Mo*+ the anisotropy field at the absolute zero for the 
odd isotopes is about 160 gauss larger than it is for the even isotopes. This 
will be noticed macroscopically by a change in the antiferromagnetic resonance 
frequency, in the field at which spin flipping takes place, in the Néel tem- 
perature, and in other macroscopic phenomena into which H, enters. 

The case of Sm** is much more extreme. Here the spin-Hamiltonian is as 
for Mot, but with g~ 0.6, A ~ 0.006 cm—', B ~ 0.025 cm, I = 7/2, or 
0. The even isotopes should behave rather like an ordinary isotropic anti- 
ferromagnetic as before. The interaction field can be estimated to be about 
30 gauss (using JJ ~ yu/d*); this produces a splitting of about 8X10-*cm™, 
which is much less than the hyperfine splitting. Thus for the odd isotopes of 
samarium the hyperfine interaction is the most important feature of the anti- 
ferromagnetic state. It is not very convenient here to talk about an anisotropy 
field (which in this case would be very much larger than the interaction 
field). Since the splitting in the interaction field is only about 8X10-' cm", 
a Néel temperature of not more than 1 or 2X10-°°K is to be expected. At 
such temperatures, only the lowest pair of hyperfine levels is appreciably 
populated, and other hyperfine levels should be ignored in treating the anti- 
ferromagnetic state, in a manner analogous to the treatment of crystalline 
field splittings in more familiar substances. If A were greater than B, it 
would be appropriate to use an “Ising model” to treat this problem. 

These two examples, Mo*+ and Sm**+, are not really representative of 
magnetic substances, since it is not usual to find an isotropic g-factor com- 
bined with an anisotropic hyperfine coupling. Usually both the g-factor and 
the hyperfine coupling are anisotropic. In this case, too, there is a significant 








H, = seo ~ 160 gauss. 
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difference between the antiferromagnetic properties of various isotopes; the 
hyperfine coupling is large enough to produce such an effect in the rare earths 


and in the actinides. 
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A NOTE ON THE SCATTERING OF A PLANE ELECTROMAGNETIC WAVE 
BY A SMALL, THIN-WALLED, CYLINDRICAL DIELECTRIC TUBE* 


J. Y. Wonc 


The problem of scattering of a plane electromagnetic wave by a hollow 
cylindrical dielectric tube, as described in this note, was undertaken as a result 
of model UHF shipborne antenna studies carried out at the Radio and Elec- 
trical Engineering Division laboratories of the National Research Council for 
the Royal Canadian Navy. In such applications, the antenna is usually a 
vertical dipole located near the top of the ship’s pole mast in order to obtain 
an omnidirectional pattern in the azimuthal plane. From our studies it was 
found that in addition to the adjacent superstructure, the pole mast itself 
contributed significantly to the pattern deterioration. Measurements showed 
that a more uniform pattern could be achieved by replacing the hollow metal 
pole with one of dielectric. 

In this note, the amplitude coefficient of the scattered field is derived for 
the case of a plane electromagnetic wave normally incident on a small thin- 
walled cylindrical dielectric tube. The results are a special case of the general 
problem of scattering by a coaxial cylinder of arbitrary constants which has 
been treated by Tang (1956). Figure 1 illustrates the geometrical arrangement 
used in the analysis. Following Tang’s notation, the scattered field in region (1) 
is of the form: 


co 


(1) Et = >> i" A, Hy” (kip) exp (ind) 


where k? = wye+iwop and A, is the scattering coefficient. 

Let ky = k3 = k, and in region (2), o2 = 0. The resulting configuration 
becomes a lossless dielectric tube having inner and outer radii of a and b 
respectively. Furthermore, it is assumed that @ and b are small in terms of 
wavelength and that a/b > 0.9. Introducing these restrictions in the general 
expression for the scattering coefficient (Eq. 22, Tang), and employing the 


*Issued as N.R.C. No. 4981. 
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INC 





Tyree y 


Fic. 1. Coaxial cylinder with arbitrary media constants. 


small argument approximation for the cylindrical functions, the following 
result is obtained for the zero-order term of the scattering coefficient, 


_ _ix,2,S(2_,)(b_2)\ 
(2) eS (e-1)(3-$) 


= —ER* (b°—a") (@-1). 


The other coefficients are of higher powers in kb and hence for this case are 
vanishingly small. As to be expected, the result given in Eq. (2) implies that the 
small dielectric tube acts as a line source, i.e., a uniform scatterer. 

It is interesting to compare our result with that obtained by Wait (1955) 
for a small solid dielectric cylinder. The zero-order term of the scattering 
coefficient for a solid cylinder is given by 

in 


, Pa 2) 


(3) ay = 
where ¢ is the dielectric constant of the cylinder material. By equating eqs. 
(2) and (3) one obtains an expression for the ‘‘equivalent radius” of the tube 
in terms of a solid dielectric cylinder. Note that the difference in signs for 
the two coefficients occurs because of the difference in the positive direction 
of the azimuthal angle of incidence. There results an expression for the 
equivalent radius, 


(4) Qequ = /(b°—a’). 
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LOW ENERGY PROTON FLUX AT SEA LEVEL* 


I. B. McDIarRMIp 


In the course of an experiment carried out to search for a 500 m, particle, 
the vertical proton flux in the momentum interval 716-828 Mev/c, at latitude 
57° N. and average atmospheric pressure 1002 mb has been measured. Such 
a measurement was thought worthwhile since considerable spread exists in 
previous measurements of the flux around 700 Mev/c. (For a summary of the 
work up to 1955 see Ogilvie 1955.) The present experiment was relatively free 
from uncertainties in the identification of protons, and had the advantage of 
a relatively small acceptance angle (half angle in one plane of 9.8° and 5.8° 
in a perpendicular plane). 

The experimental arrangement shown in Fig. 1 is similar to one described 
previously (McDiarmid 1958). The cloud chamber contained two plastic 
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Fic. 1. Experimental arrangement. 


scintillators and nine Al plates each 2.16 g/cm? thick. The chamber was 
expanded and the scintillator pulse heights recorded when the following 
selection criterion was satisfied: 


A+C+(T > 12 v)+(B > l4v)-D 


where 7’ and B refer to the pulse heights from top and bottom scintillators 
respectively. (The most probable pulse size due to a relativistic ~-meson 
travelling vertically was 8 v, and that due to a proton stopping in plate 3 
was 23 v and 28 v in top and bottom scintillators respectively.) This selects 


*Issued as N.R.C. No. 4989. 
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protons stopping in one of the plates in the chamber with an efficiency of very 
nearly 100%, but u-mesons stopping are detected with an efficiency of only 
6%. The apparatus operated 1500 hours from March to June 1958 during which 
time 270 particles were observed to stop in the bottom scintillator and 504 in 
the Al plates. The tracks of stopping particles were reprojected and the pulse 
heights corrected for variation of scintillator response with position as des- 
cribed previously (McDiarmid 1958). The corrected pulse heights and the 
residual ranges were used to identify protons as shown in Fig. 2 in which the 
top scintillator pulse height is plotted against the bottom scintillator pulse 
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Fic. 2. Pulse height divided by path length in top scintillator vs. corresponding quantity 
in bottom scintillator for particles stopping in plate 3. 


height for particles stopping in plate 3. In this way 100 particles stopping in 
the bottom scintillator and 206 particles stopping in the nine Al plates were 
identified as protons. A less satisfactory separation of protons and u-mesons 
was also made from a visual estimate of the ionization of the stopping particle 
in the last few compartments. 

For protons stopping in the lower scintillator the selection system required 
that on entering this scintillator the particles had, on the average, an energy 
> 20.5 Mev to be detected. This energy corresponds to a residual range in the 
scintillator of 0.44 g/cm? and this thickness had to be added to the total 
amount of absorber above the bottom scintillator (this included the top 
scintillator, the chamber walls, counter trays, supports, and roof) to obtain 
the minimum range of the protons detected. The range limits in g/cm? Al 
equivalent* are given in Table I for protons stopping in the bottom scintillator 
and in the Al plates. Also given in Table I are the corresponding momentum 
limits and the absolute proton intensities obtained using a solid angle X area 
X time factor, QAt = 2.55107 cm? sec, which was corrected for counter wall 


*The range energy tables in U.C.R.L. Report 2301 were used to convert thicknesses of the 
various materials to equivalent thicknesses of Al. 
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thickness and chamber dead time. A mean absorption length of 170 g/cm? 
was used to correct for proton absorption in the material above positions in 
the chamber where particles were observed. No scattering corrections have 


been applied since these were estimated to be less than 7%. 


TABLE I 





Number of Range interval, Momentum interval, Intensity particles/ 
protons observed g/cm? Al equiv. Mev/c cm?/sec/sterad/(Mev/c) 
100 45 .5-50 716-739 2.4+.3X1077 
206 50 -69.4 739-828 1.3+.1X1077 


The intensities obtained here are plotted in Fig. 3 along with those given 
by other authors in the same momentum region. The point at 783 Mev/c is in 
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Fic. 3. Proton intensity vs. momentum. 1, Ballam and Lichtenstein (1954); 2, Merkle, 
Goldwasser, and Brode (1950); 3, Mylroi and Wilson (1951); 4, York (1952) ; 5, Filthiith (1955) ; 
6, Ogilvie (1955). 


agreement with values reported by Mylroi and Wilson (1951), York (1952), 
and Ballam and Lichtenstein (1954) while the point at 727 Mev/c supports 
the higher intensity values, being at least a factor of two higher than the 
measurements of Filthiith (1955). 


The author wishes to thank Dr. D. C. Rose for the interest he has taken 
in this work. 
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LETTERS TO THE EDITOR 





Under this heading brief reports of important discoveries in physics may be published. These 
reports should not exceed 600 words and, for any issue, should be submitted not later than six weeks 
previous to the first day of the month of issue. No proof will be sent to the authors. 


Nuclear Orientation of Cobalt-60 in Antiferromagnetic Co(NH,)2(SO,)2.6H,O 


We have prepared a single crystal of cobalt ammonium sulphate, containing about 5yc 
of Co-60, and weighing about 0.25 g. This crystal was cooled to approximately 0.05° K, in 
contact with potassium chromium alum, by adiabatic demagnetization. y-Ray counters were 
placed along the AK; and Ae axes of the crystal to count the y-rays emitted by the Co-60, 
and to determine the anisotropy of emission of the y-rays. (For a detailed description of the 
crystallography of this salt, and of the techniques involved in nuclear orientation experiments, 
see Bleaney et al. (1954).) Anisotropies of up to 5}% were measured. A very similar experiment 
to this, using the antiferromagnetism of cobalt ammonium sulphate to produce nuclear 
alignment, was tried by a group at Leiden, Holland, in 1951, du: ing the time when workers 
in this field were searching for a method of achieving nuclear alignment. This group reported 
a negative result (Poppema 1954); however, Poppema does not give enough detail for us 
to try to explain why a negative result should have been obtained. In view of this, we believe 
it is of value to describe our experimental arrangement and results. 
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Fic. 1. After 180 minutes, the sample was warmed up to 1.3° K to obtain a normalizing count, which is shown, 
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‘To one side of a copper disk, | cm diameter, were soft-soldered 100 No. 36 B&S copper 
wircs; surface area about 70 cm?. This disk was placed in a pill press, potassium chromium 
alum powder was poured on to the wires, and the whole was compressed. This formed a cylinder 
of chrome alum, 1 cm diameter X 4 cm long, with a copper disk on one end. The wires, which 
stretched the whole length of the cylinder, served to conduct heat between the copper disk 
and the remote corners of the chrome alum. The crystal of cobalt ammonium sulphate was 
sandwiched between the copper disks of two such chrome alum cylinders, and thermal 
contact was provided between the disks and the crystal by means of Apiezon grease N. The 
whole sample was cooled by adiabatic demagnetization from a field of 21,000 oersted and an 
initial temperature of 1.3° K. The field was applied along the a-axis of the cobalt ammonium 
sulphate crystal, and was turned off slowly in a period of 4 minutes to avoid eddy current 
heating in the copper. The anisotropy of the y-rays was measured in the usual way (Bleaney 
et al. 1954), and the temperature of the chrome alum was measured by the ballistic method. 

The results of one experiment are shown in Fig. 1. It is seen that the anisotropy grows 
from zero, or almost zero, to about 5.5%, and then decreases as the whole sample warms up. 
If we assume that there is thermal equilibrium during the period while the anisotropy is decreas- 
ing (a debatable assumption which we do not wish to justify in this short communication) 
we can obtain values of anisotropy as a function of temperature. Such values are plotted in 
Fig. 2: the solid curve is the value of the anisotropy obtained from a paramagnetic Tutton 
salt (87% Zn, 12% Cu, 1% Co) Rb2(SO,)2.6H2O (Bleaney et al. 1954). It is interesting to note 
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that the anisotropy is very nearly the same as is obtained in a paramagnetic salt. We can 
conceive of this if we admit that the electron “spin"’ of each ion is constrained to point along 
a definite crystal axis, the tetragonal axis of the crystalline field (Bleaney et a/. 1954). The 
crystalline field and hyperfine coupling produce such a large anisotropy energy that the spin 
may point in either direction along this axis, but may not point at an angle to this axis. In 
the paramagnetic state, the direction of the spin of any one ion is decided at random, while 
in the antiferromagnetic case the direction is determined by the interaction. The anisotropy 
of y-radiation does not depend on the direction of the electron spin. 

Cobalt ammonium sulphate has been studied by Garrett (1951). According to his work, the 
crystal should have cooled to about 0.05° K immediately after demagnetization, and an 
anisotropy of 5.5% should have been noticed from the start. At the present time, we can offer 
no explanation why this is not so. 


We wish to acknowledge the support of the National Research Council of Canada for this 
work. One of us (M.A.R.L.) is grateful for the award of an N.R.C. studentship during the 
tenure of which this work was carried out. 
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Announcement on the International Yard and Pound 


It would be appreciated if you would publish the enclosed announcement on the international 
yard and the international pound. The text is self-explanatory but it may be of interest to 
recall the following facts of local concern to Canada. 

In 1951 ‘‘An Act Respecting Units of Length and Mass”’ (15 George VI, Chapter 31) defined 
the Canadian yard as 0.9144 meter, which is in exact agreement with the joint announcement. 
The same Act defined the Canadian pound as 0.45359243 kilogram which is slightly discrepant 
with the ratio being established for the International pound. This difference is insignificant 
in all measurements except those of the very highest precision. Nevertheless steps are being 
taken toward the amendment of the Canadian legislation to make the Canadian pound 
0.45359237 kilogram and so legally equal to the International pound. The ratio’ 0.9144 was 
chosen because it leads to a yard which is approximately half-way between the present Imperial 
(United Kingdom) yard and the U.S. yard. It has the added advantage of making one inch 
exactly equal to 2.54 centimeters. The ratio used in the Canadian legislation of 1951 to relate 
the Canadian pound to the International kilogram was the originally accepted ratio between 
the pound and the kilogram. Through the years the Imperial (United Kingdom) pound and 
the U.S. pound have become very slightly different and the ratio used in the joint announce- 
ment makes the International pound approximately the average between the Imperial pound 
and the U.S. pound. 


ANNOUNCEMENT 
The Directors of the following standards laboratories: 


Applied Physics Division, National Research Council, Ottawa, Canada 
Dominion Physical Laboratory, Lower Hutt, New Zealand 

National Bureau of Standards, Washington, United States 

National Physical Laboratory, Teddington, United Kingdom 

National Physical Research Laboratory, Pretoria, South Africa 
National Standards Laboratory, Sydney, Australia 


have discussed the existing differences between the values assigned to the yard and to the 
pound in different countries. To secure identical values for each of these units in precise 
measurements for science and technology, it has been agreed to adopt an international yard 
and an international pound having the following definitions: ’ 


the international yard equals 0.9144 metre; 
the international pound equals 0.453 592 27 kilogramme. 


It has also been agreed that, unless otherwise required, all non-metric calibrations carried 
out by the above laboratories for science and technology on and after July 1, 1959, will be 
made in terms of the international units as defined above or their multiples or submultiples. 


RECEIVED DECEMBER 16, 1958. L. E. Howvet1 
DIVISION OF APPLIED PHysICcs, 

NATIONAL RESEARCH COUNCIL, 

Ottawa, CANADA. 
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points, periods, and stippled dots must be solid black circles large enough to be reduced if 
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